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ABSTRACT 

We present a set of cluster models that link the present-day properties of clusters to the processes that 
govern galaxy formation. These models treat the entropy distribution of the intracluster medium as its 
most fundamental property. Because convection strives to establish an entropy gradient that rises with 
radius, the observable properties of a relaxed cluster depend entirely on its dark-matter potential and the 
entropy distribution of its uncondensed gas. Guided by simulations, we compute the intracluster entropy 
distribution that arises in the absence of radiative cooling and supernova heating by assuming that the 
gas-density distribution would be identical to that of the dark matter. The lowest-entropy gas would then 
fall below a critical entropy threshold at which the cooling time equals a Hubble time. Radiative cooling 
and whatever feedback is associated with it must modify the entropy of that low-entropy gas, changing 
the overall entropy distribution function and thereby altering the observable properties of the cluster. 
Using some phenomenological prescriptions for entropy modification based on the existence of this cooling 
threshold, we construct a remarkably realistic set of cluster models. The surface-brightness profiles, mass- 
temperature relation, and luminosity-temperature relation of observed clusters all naturally emerge from 
these models. By introducing a single adjustable parameter related to the amount of intracluster gas 
that can cool within a Hubble time, we can also reproduce the observed temperature gradients of clusters 
and the deviations of cooling-flow clusters from the standard luminosity-temperature relation. 

Subject headings: cosmology: theory — galaxies: clusters: general — galaxies: evolution — 
intergalactic medium — X-rays: galaxies: clusters 



1. INTRODUCTION 

The hot gas in clusters of galaxies seems like it should 
be easy to understand. Because of the relatively low ratio 
of baryons to dark matter, the potential well of a clus- 
ter should be dark-matter dominated. The dynamical 
time within a cluster potential is shorter than a Hubble 
time, so most clusters should be relaxed. Also, the cool- 
ing time of the vast majority of intracluster gas is longer 
than a Hubble time. It would appear that cluster struc- 
ture ought to be scale-free, as long as the shape of a clus- 
ter's potential well does not depend systematically on its 
mass. If that were the case, then the global properties 
of clusters, such as halo mass, luminosity-weighted tem- 
perature, and X-ray luminosity, would scale self-similarly 
(Kaiser 1986). In particular, the gas temperature would 
scale with cluster mass as T cx M 2 ^ 3 and the bolometric 
X-ray luminosity would scale with temperature as L oc T 2 
in the bremsstrahlung-dominated regime above ~ 2 keV. 
Indeed, numerical simulations that include gas dynamics 
but exclude non-gravitational processes such as radiative 
cooling and supernova heating produce clusters that obey 
these scaling laws (e.g., Evrard, Metzler, & Navarro 1996; 
Bryan & Norman 1998; Thomas et al. 2001b). 

Real clusters are not so simple. We have known for a 
decade that the observed luminosity-temperature relation 
is closer to L oc T 3 (e.g., Edge & Stewart 1991), indi- 
cating that non-gravitational processes must influence the 
density structure of a cluster's core, where most of the lu- 
minosity is generated (Kaiser 1991; Evrard & Henry 1991; 



Navarro, Frenk & White 1995; Bryan & Norman 1998). 
The mass-temperature relation, on the other hand, seemed 
like it ought to be more fundamental and less sensitive to 
non-gravitational effects. Yet, observations collected over 
the last few years indicate that this relation also disagrees 
with both the scale-free predictions and simulations that 
exclude non-gravitational processes (Horner, Mushotzky, 
& Scharf 1999; Nevalainen, Markevitch, & Forman 2000; 
Finoguenov, Reiprich, & Bohringer 2001; Xu, Jin, & Wu 
2001). These results derive mostly from resolved X-ray 
and temperature profiles coupled with the assumption of 
hydrostatic equilibrium, but they do seem consistent with 
gravitational lensing measurements (Allen, Schmidt, & 
Fabian 2001). 

Understanding the scaling properties of clusters is of 
broad importance because these scaling laws are integral 
to determination of cosmological parameters. For exam- 
ple, the number density of clusters in the local universe 
reflects the amplitude of matter perturbations on ~ 20 
Mpc scales (e.g., Henry & Arnaud 1991; White, Efs- 
tathiou, & Frenk 1993; Viana & Liddle 1996; Eke, Cole k 
Frenk 1996; Kitayama & Suto 1997; Oukbir & Blanchard 
1997; Henry 2000). In order to measure this amplitude, 
one needs accurately measured cluster masses. Gravita- 
tional lensing measurements of cluster mass have grown 
rapidly in sophistication during the last few years, but X- 
ray determined cluster temperatures, converted to cluster 
masses using the mass-temperature relation, remain the 
most commonly used mass measurement for cosmological 
studies. Thus, any inaccuracies in the mass-temperature 
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relation propagate into uncertainties in cosmological pa- 
rameters derived from clusters (e.g., Voit 2000). A num- 
ber of authors have recently re-evaluated the perturba- 
tion amplitude implied by cluster observations using a va- 
riety of techniques that circumvent the standard mass- 
temperature relation (Rciprich & Bohringer 2001; Van 
Waerbeke et al. 2001; Seljak 2002; Viana, Nichol, & Lid- 
die 2002). The results are disconcertingly different from 
those based on the standard relation, generally implying a 
much lower power-spectrum amplitude. Clearly a better 
understanding of clusters is needed, and ideally, we would 
like to know how the present-day properties of clusters are 
linked to the physics of galaxy formation. 

The physical processes most likely to break the expected 
self-similarity of clusters are heating by supernovae or ac- 
tive galactic nuclei and radiative cooling, although other 
possibilities such as magnetic pressure or cosmic-ray pres- 
sure have not been ruled out. A large number of stud- 
ies have investigated how the structure of the intracluster 
medium changes when these processes are included. Be- 
cause treating all of these processes in sufficient detail is 
extremely complicated, various approximations have been 
employed. These include numerical simulations that ad- 
minister a particular amount of energy input at early times 
or adopt some other phenomenological prescription for 
injecting energy (e.g., Metzler & Evrard 1994; Navarro, 
Frenk, & White 1995; Metzler & Evrard 1997; Bialek, 
Evrard, & Mohr 2001; Borgani et al. 2001a; da Silva et 
al. 2001; Bryan & Voit 2001), analytical models for spher- 
ically symmetric accretion of preheated gas (e.g., Cava- 
liere, Menci & Tozzi 1999; Tozzi & Norman 2001), and 
semi-analytic models that consider the merger and star 
formation history of clusters (e.g., Wu, Fabian, & Nulscn 
1998, 2000; Balogh, Babul, & Patton 1999; Bower et al. 
2001; Babul, Balogh, Lewis, & Poole 2002). Most of these 
papers have focused on heating of the gas, but recent work 
has also examined the possibility that radiative cooling is 
primarily responsible (e.g., Knight & Ponman 1997; Sug- 
inohara & Ostriker 1998; Pearce et al. 2000; Bryan 2000; 
Muanwong et al. 2001; Thomas et al. 2001a; Lewis et al. 
2001). Two broad conclusions emerge from such studies: 
(1) radiative cooling without feedback probably locks too 
high a fraction of baryons into condensed objects, and (2) 
if heating is uniformly distributed, the amount required to 
produce realistic clusters is quite large, probably exceed- 
ing ~ 1 keV per baryon. Introducing this much energy 
through supernova heating strains the bounds of plausibil- 
ity, given how many metals are observed in the intracluster 
medium (e.g., Balogh, Babul, & Patton 1999; Valageas & 
Silk 1999; Kravtsov & Yepes 2000; Bower et al. 2000; 
although see Loewenstein 2000). 

Here, we revisit an approach pioneered by Bower (1997), 
focusing on entropy as the most fundamental character- 
istic of the intracluster medium. Analyzing clusters in 
terms of entropy has two advantages. The first is that 
convection within clusters acts as an entropy-sorting de- 
vice, shuttling low-entropy material to the cluster core and 
high-entropy material to the cluster's outskirts. Convec- 
tive stability is achieved when specific entropy becomes 
a monotonically increasing function of radius. Thus, if 
we can determine the entropy distribution of the gas in 
a cluster and we know the structure of its dark-matter 



halo, we can predict the structure of that cluster in its 
relaxed state. In the idealized cases that we will consider, 
we will assume that clusters are spherically symmetric, 
dark-matter dominated, in hydrostatic equilibrium, and 
in convective equilibrium. The second advantageous prop- 
erty of entropy is that, in the absence of non-gravitational 
processes, the entropy distribution of a cluster can be read- 
ily computed using either hydrodynamical simulations or 
spherically symmetric analytical models. Then, the effects 
of non-gravitational heating and radiative cooling can be 
implemented as modifications of that entropy distribution. 

In Voit & Bryan (2001), we introduced a simple scheme 
for computing cluster models, taking advantage of these 
properties of the intracluster entropy distribution to show 
how heating and cooling conspire to produce a cluster L-T 
relation quite similar to what is observed. This paper ex- 
plores those modified-entropy models in more detail, show- 
ing that the observable properties of these model clusters 
are quite similar to those of real clusters on many counts. 
Section 2 describes how we construct those models and an- 
alyzes how cluster properties depend on both halo concen- 
tration and entropy modification. Section 3 uses relations 
between halo mass and halo concentration to define sets 
of cluster models determined entirely by the underlying 
cosmology. We then compare the model clusters with ob- 
servations of the surface-brightness profiles of clusters, the 
temperature gradients of clusters, the mass-temperature 
relation, and the luminosity-temperature relation. In ev- 
ery case, our models reproduce the observations, with lit- 
tle need for adjustable parameters. In section 4, we ex- 
plore how the present-day intracluster entropy distribution 
emerges from hierarchical cosmological models. Section 5 
summarizes our conclusions. 

2. MODIFIED ENTROPY MODELS 

In this section we develop a family of modified-entropy 
models of clusters based on two parameters, one speci- 
fying the concentration of the cluster's dark-matter halo 
and the other specifying the entropy level at which cool- 
ing and feedback modify the intracluster entropy distri- 
bution. We begin by casting the equilibrium equations 
in dimcnsionlcss form and selecting appropriate boundary 
conditions. Then we assume that the unmodified entropy 
distribution determined by gravitational processes alone 
would produce a gas density distribution similar to the 
dark-matter density distribution. Because the innermost 
gas in such a cluster can cool in much less than a Hubble 
time, we proceed to investigate several potential modifica- 
tions of the entropy distribution by radiative cooling and 
subsequent feedback. In each case, we show how both en- 
tropy modification and the shape of the underlying dark- 
matter halo affect the surface-brightness profile of a clus- 
ter, its luminosity-weighted temperature, its temperature 
gradient, and its X-ray luminosity. Two important results 
emerge from this analysis: (1) the temperature of a cluster 
of a given mass is determined primarily by halo concen- 
tration and is affected only modestly by entropy modifica- 
tion, and (2) altering the intracluster entropy distribution 
at the scale set by radiative cooling leads to an Lx-T rela- 
tion that has the observed slope at z = and that evolves 
very little with time. 
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2.1. Equilibrium Structure 

The equilibrium structure of the intracluster medium 
is determined by three things, the gravitational potential 
of the cluster's dark matter, the entropy distribution of 
the cluster's gas, and the confining pressure at the outer 
boundary of the cluster. In order to solve for that struc- 
ture, one must integrate the equations of hydrostatic equi- 
librium and gas mass conservation, 
dP GM(< r) 

15 P I 1 ) 



dr 



dM g 
dr 



Airr p , 



(2) 



using the equation of state P — Kp 1 , where 7 is the adi- 
abatic index and K specifies the adiabat of the gas. In 
this system of equations, M(< r) is the total mass within 
radius r, M g is the gas mass within that radius, and the 
other symbols have their usual meanings. Throughout this 
paper, we will assume that the gas mass is gravitationally 
negligible and that 7 = 5/3, as appropriate for an ideal 
monatomic gas. 

In a sufficiently relaxed cluster, convection ensures that 
K(r) monotonically increases with radius. Thus, if one 
knows M g (K), the entropy distribution of the intraclus- 
ter medium expressed in terms of the mass of gas with 
Pp-5/3 < K} then one can use the inverse relation K(M g ) 
to solve the equilibrium equations that determine the clus- 
ter's structure. Because this crucial quantity K is so 
closely related to specific entropy, we will often refer to it 
as the "entropy" of the gas, even though the standard ther- 
modynamic entropy per particle for an ideal monatomic 
gas is s + const. 

2.1.1. Dimensionless Form 

We can gain insight into the key parameters that control 
the structure of the intracluster medium by investigating a 
family of dimensionless models for clusters in hydrostatic 
and convective equilibrium. Because the virial radius of 
a cluster lies close to the radius r 2 oo within which the 
mean matter density of the cluster is 200 times the crit- 
ical density p cr , we elect to represent radii in terms of 
f = r/r 2 oo and gas density in terms of p = p/fbP200, where 
P200 = 200p cr and = 0.02/fl m h 2 is the fractional con- 
tribution of baryons to the total mass of the universe. One 
can then rewrite the equilibrium equations in dimension- 
less form as follows: 



dP „M„ 
dr r z 

dr 



3f 2 p , 



(3) 
(4) 



where M200 = (47r/3)r| 00 p200 is the virial mass, T 2 oo = 
GM 2 ooM TO p/2r , 20o is the temperature in energy units 1 

of the corresponding singular isothermal sphere, P = 
P/[T2oofbP20o(^n l p)^ 1 ] is the dimensionless pressure, M = 
M(< r)/M 2 oo is the dimensionless mass within f, and 
f g = M s /(/;,M2oo) is the fraction of a cluster's baryons 
in the intracluster medium within radius f. The corre- 
sponding dimensionless entropy distribution is K(f g ) = 
K /[T2oo(.fbP2oa)~ 2 ^ 3 (p m p)~ 1 ], and we will also make use 
of the quantity T = P p^ 1 puipT^- 

1 Boltzmann's constant k is absorbed into T throughout the paper. 



General solutions to these equations can be found for 
a few idealized cases. For example, if the potential is a 
singular isothermal sphere and the gas is also isothermal, 
then we have P oc p oc f~ 2 , f g oc f, and K oc f 4 / 3 for a 
monatomic ideal gas. If instead the mass profile follows 
the form of Navarro, Frenk & White (1997; NFW here- 
after), with M oc [ln(l + cf) — cf(l + cf) _1 ] where c is the 
concentration parameter, then the equation of hydrostatic 
equilibrium becomes (see Wu, Fabian, & Nulscn 2000) 



dP 

df 



= 2p 



ln(l + c) - 



- 1 d 


"ln(l + cr)~ 


dr 


f 



(5) 



Given an isentropic gas in which Pp~ 5 ^ 3 = K = const, 
we thus obtain p oc [f _1 ln(l + cf) + C] 3 / 2 , where C is a 
constant of integration. 

2.1.2. Boundary Conditions 

A particular solution to the equilibrium equations can 
be found by choosing P and f g at the origin and inte- 
grating outwards. The gas mass fraction at the origin is 
always zero by definition, but there is some freedom in the 
choice of P = P(0). Any given choice of Pq corresponds 
to a unique ICM structure with a unique luminosity and 
temperature profile, but which solutions are the physical 
ones? 

Here is where the confining pressure comes into play. 
It is generally assumed that accreting matter confines the 
ICM in the neighborhood of r 2 oo- The pressure at r 2 oo 
would then be determined by the ram pressure of the in- 
falling matter. Now suppose that K(M g ) has been modi- 
fied by non-gravitational processes. For example, heating 
might drive gas out of the cluster potential, pushing the ac- 
cretion shock to a radius beyond r 2 oo- Alternatively, cool- 
ing and condensation might reduce the total gas mass of 
the intracluster medium so that the accretion shock moves 
inward. In either case the radius at which the integration 
of hydrostatic equilibrium should terminate differs from 
r 2 oo- 



So what is the termination radius 



at which 



/ g /bM 2 oo equals the total uncondensed gas mass, and what 
is the pressure at that radius? In a freely falling accretion 
flow with a constant mass flux, the density should vary 
as vj^r^ 2 , where Vff oc f _1 / 2 '\/ln(l + cf) is the free-fall 
velocity from f > 1 in an NFW potential. We therefore 
assume that the accretion pressure scales with radius as 
P a cc(f) oc f _5 / 2 ^/ln(l + cf ), and we normalize that pres- 
sure to equal the unmodified NFW value (see § 2.2) at 
f = 1. Various choices of Pq will lead to pressure profiles 
that intersect -P QCC (f) at different values of f max , each cor- 
responding to a unique value of fg(f max ) determined by 
a particular integration. The correct solution is the one 
that gives the proper value for the uncondensed gas mass 
at the radius where the gas pressure equals the accretion 
pressure. 

In practice, the physically reasonable solutions are not 
particularly sensitive to the outer boundary condition. In 
the models that follow, applying the boundary condition 
P = Pace at f = 1 would lead to similar results. How- 
ever, the value of f g derived at f = 1 using this boundary 
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condition is not physically significant and can even be in- 
consistent with models in which cooling and condensation 
significantly reduce the maximum value of f g . One exam- 
ple of a physically unreasonable boundary condition is re- 
quiring f g = 1 at f = 1 . Such a requirement does not allow 
intracluster gas to expand beyond r 2 oo even when heating 
substantially raises the intracluster entropy. In the limit of 
extreme heating, this condition therefore artificially boosts 
the central pressure, density, and temperature, leading to 
an unphysically large cluster luminosity. 

2.1.3. Integrated Characteristics 

Once a correct solution has been identified, one can com- 
pute two quantities that correspond to the dimensionlcss 
luminosity and luminosity-weighted temperature: 

ff max 

L = 3 p 2 f 2 df (6) 

Jo 

T] um - 3L- 1 / fp 2 f 2 df . (7) 
Jo 

As long as the cooling function A(T) remains sufficiently 
constant within the gas contributing the bulk of the lu- 
minosity, the actual luminosity-weighted temperature is 
Ti um « X^oo^lum, and the bolometric X-ray luminosity is 

L " T r ™ (5) ^ A(Tlum)i ' (8) 

where n e = {n e /p)f b p 2 oo = 1-2 x 1CT 4 (^m/0.33)" 1 cm~ 3 . 
Scaling to a 10 keV cluster at z = gives 

^ (9 .0x 1 0^-»e Ig ,-,(^) S,2 (5|)" 2 

x A 23 (T lum )L , (9) 

where A 23 (Ti um ) = A(21 lim )/(l()- 23 ergcm 3 s- 1 ). 

The emission-measure profile and emissivity-wcightcd 
temperature profile are also of interest. Thus, we define 
the following dimensionlcss analogs as functions of the pro- 
jected radius f±: 

A f"\f ^ max ^ ±_ 

S(f±) = 2 / p 2 dl (10) 

Jo 

T ± (f ± ) = 2S~ 1 / Tp 2 dl , (11) 

Jo 

where I = \/r — fj_ is the dimensionless distance along the 
line of sight. 

2.2. Unmodified Solutions 

Before exploring the consequences of modifying the en- 
tropy of intracluster gas, we need to know what the en- 
tropy distribution would be if it were not modified. To 
simplify matters, we will assume that the underlying dark 
matter profile of the cluster is of NFW form with some con- 
centration parameter c. This parameter typically ranges 
from c ~ 5 for hot clusters toc~ 10 for groups of galaxies 
(see § 3.1). If the intracluster medium were collisionless, 
then its density profile would be identical to that of the 
dark matter. Thus, we define the unmodified entropy dis- 
tribution K (M g ) of a cluster of concentration c to be that 
of gas in hydrostatic equilibrium in the cluster potential 
with a density profile identical to that of the dark matter. 
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Fig. 1. — Intracluster density profiles in dimensionless units. The 
dotted line shows the matter density (p) of an NFW profile of con- 
centration c = 8 in units of the mean density within 7-200 as a 
function of radius (f) in units of r200- The dashed line shows the 
density profile of a simulated cluster whose dark matter halo is well 
approximated by an NFW halo with c = 8. Outside of f = 0.1, 
these curves closely correspond. The solid lines show how the run 
of density with radius changes as the entropy distribution of an un- 
modified c = 8 NFW halo is truncated at progressively larger values 
of the dimensionlcss entropy K c . For these cases, p is the gas density 
in units of the mean density of the unmodified gas. Notice that the 
model with K c = 0.1 is quite similar to the gas density distribution 
from the simulation. 

Numerical experiments reveal that the entropy distri- 
bution of intracluster gas in clusters simulated without 
radiative cooling or supernova heating is quite similar to 
the NFW form throughout most of the cluster but is el- 
evated above NFW within the cluster's core. Figure 1 
shows the density distribution of a cluster simulated with 
an adaptive-mesh refinement code (Norman & Bryan 1998; 
Bryan 1999) for the Santa Barbara cluster comparision 
project (Frenk et al. 1999). The dotted line labeled 
"NFW" shows an NFW density distribution with c = 8, 
which is very close to the underlying dark-matter distri- 
bution of the cluster. The dashed line labeled "simula- 
tion" shows the gas-density distribution from that same 
simulation. Agreement is quite close outside 10% of the 
virial radius, but the gas density levels off within that ra- 
dius, implying that the lowest-entropy gas has a somewhat 
higher entropy than assumed in our unmodified distribu- 
tion Kq. However, this discrepancy at the low-entropy 
end of the distribution is inconsequential to the modificd- 
entropy models that follow because the entropy of the dis- 
crepant gas will always be subject to further modification. 

The unmodified distribution is also somewhat dis- 
crepant with simulations at the high-entropy end. Assum- 
ing hydrostatic equilibrium and an NFW gas-density law 
near the virial radius leads to gas temperatures that are 
20-30% higher than those in our fiducial simulation. Thus, 
our unmodified entropy levels at large radii are slightly 
higher than what a numerical simulation would produce. 
However, in order to properly reproduce the density pro- 
file near the virial radius in a hydrostatic model, we must 
retain these elevated entropy levels. This overestimate at 
the cluster's outskirts has virtually no effect on emissivity- 
weighted global quantities, such as L and Ti um , but if one 
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is interested in the temperature of cluster gas near r2oo, 
then departures from hydrostatic equilibrium, which are 
outside the scope of our models, must be taken into ac- 
count. 

2.3. Entropy Modification 

Modification of this baseline entropy distribution is in- 
evitable because clusters simulated without radiative cool- 
ing are physically inconsistent. In general, the gas at the 
center of such a simulated cluster can radiate many times 
its thermal energy within a Hubble time, resulting in clus- 
ter luminosities that vastly exceed those observed (e.g., 
Muanwong et al. 2001; Bryan & Voit 2001). Condensa- 
tion and removal of the lowest-entropy gas from the intr- 
acluster medium must happen at some level, because that 
is how the cluster's galaxies form stars. Furthermore, the 
most massive of these stars must explode, resulting in su- 
pernova feedback. Condensed gas accreting onto an active 
galactic nucleus can provide additional feedback (Heinz, 
Reynolds, & Begelman 1998; Kaiser & Alexander 1999; 
Quilis, Bower, & Balogh 2001; Reynolds, Heinz, & Begel- 
man 2002; Bohringer et al. 2002). All of these processes 
will modify the entropy distribution of the intracluster 
medium, but modeling them in detail is a daunting task. 

Because the physics of feedback is so complex, we shall 
adopt a highly simplistic phcnomenological approach to 
entropy modification. Instead of trying to model all the 
consequences of heating and cooling we will restrict our 
investigation to three qualitatively different modifications 
of the entropy distribution K : (1) truncation of the dis- 
tribution, (2) shifting of the distribution, and (3) radia- 
tive losses from the distribution. Each of these modifica- 
tions depends on a single entropy threshold parameter K c , 
which we take to be the entropy at which the cooling time 
of gas of temperature T 2 oo equals the age of the universe. 
Thus, for each type of entropy modification, we obtain a 
two-parameter family of modified-entropy models specified 
by a concentration c and a dimcnsionless entropy K c . 

As in Voit & Bryan (2001), we use an approximate cool- 
ing function for gas with a metallicity of one-third solar 
to determine the appropriate entropy threshold. Accord- 
ing to this approximation, the threshold for cooling within 
15 Gyr can be expressed as 



K r 



135 keV cm 



-2/3 



2kcvy 



C 



(12) 



/j,m p 

with C = 2/3 for T > 2kcV and C = for T < 2keV. 
Converting this threshold to dimensionless units yields 

K r . « 0.164 f lum ( 7^) ' ' . (13) 



2keV 



As long as f\ x 



1 (see § 2.5), then K c 



0.1 will be 



appropriate for the hottest clusters (<~ lOkeV) and K w 
0.5 — 1.0 will be appropriate for the coolest (<~ 0.5 keV) 
groups. 

2.3.1. Truncation 

One simple yet physically motivated way to modify the 
entropy distribution of intracluster gas is to truncate it at 
K c (Bryan 2000; Voit & Bryan 2001). Mathematically, we 



<^ 1 




Fig. 2. — Modified entropy distributions. The thick solid line 
shows how dimensionless entropy K rises with f g , the fraction of a 
cluster's baryons in the intracluster medium within radius f, for an 
unmodified c = 8 NFW halo. Gas in the unmodified halo is assumed 
to be in hydrostatic equilibrium with a density profile identical to 
that of the dark matter. The thin solid line illustrates the entropy 
distribution Kx(fg) that results when the unmodified distribution 
is truncated at K c = 0.3, corresponding to a horizontal translation 
of the unmodified distribution. The dotted line shows Ks(fg), the 
unmodified distribution shifted vertically by K c = 0.3. The long- 
and short-dashed lines show two different modifications designed to 
mimic radiative losses from the unmodified distribution (see § 2.3.3). 
These lines correspond to a reduction of K at each point of the un- 
modified distribution, followed by a horizontal translation so that 
K(0) = 0. 

define the truncated distribution to be 

K T (f g ) = K (f g + f c ) , (14) 

where f c is defined by Ko(f c ) — K c . Such a modification 
corresponds to removing all the gas with a cooling time 
less than a Hubble time and ignoring the effects of cooling 
on the rest of the intracluster medium. This type of mod- 
ification can viewed either as an extreme form of cooling, 
in which all the gas below the critical threshold condenses, 
or an extreme form of heating, in which all the gas below 
the threshold is heated to very high entropy (K 3> 1) and 
convects beyond the virial radius of the cluster. Figure 1 
shows some gas density profiles resulting from different 
values of K c . The thin solid line labeled Kt in Figure 2 
shows the truncated entropy distribution for c = 8 and 
K c = 0.3. 

Notice that the maximum value of f g in these trun- 
cated models is less than unity because some of the gas 
has been removed from the intracluster medium. The 
thick solid line in Figure 3 illustrates how the fraction 
fg.max = 1 - / c of a cluster's baryons that remain in the 
intracluster medium depends on K c . For hot clusters this 
fraction is ~90%, but for cool groups it drops below 50%. 
Nevertheless, the outer radius f max of the intracluster gas 
remains close to r 2 oo, even when f g , ma ,x drops below 0.5. 
Figure 1 shows the reason. As the entropy threshold for 
truncation rises, the central entropy of the cluster also 
rises. The resulting density profile is therefore much flat- 
ter for large values of K, enabling a smaller amount of gas 
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Fig. 3. — Maximum radius and baryon content of the intraclustcr 
medium in modified-entropy models. The thick solid line shows how 
fg,max, the fraction of a cluster's baryons residing in the intraclus- 
ter medium, depends on the entropy threshold K c in the modified 
distributions K^, and Kg3- The other lines show how r max , 

the maximum radius of the intracluster medium, depends on K c in 
each of these modified distributions. 

to fill a similar volume. Note also that the density profile 
of a modified-entropy model with c = 8 and K c = 0.1 is 
quite similar to that of a c = 8 cluster simulated without 
cooling and feedback, suggesting that the cooling thresh- 
old may have only a modest effect on the structure of the 
hottest clusters. 

2.3.2. Shift 

Another way to modify intracluster entropy is to add 
a constant term K c to the unmodified entropy distribu- 
tion (Voit & Bryan 2001). This type of modification is a 
convenient way to mimic "preheated" models for cluster 
formation that inject a fixed amount of entropy per par- 
ticle into the intergalactic medium at some early time. In 
this case, the modified entropy distribution is defined by 

K s (f g ) = k (f g ) + k c . (15) 

The dotted line in Figure 2 shows this distribution for 
c = 8 and K c — 0.3. Because entropy has been added 
but no gas has been removed, the outer radius for mod- 
els with a shifted entropy distribution extends beyond the 
usual virial radius (see Figure 3). Qualitatively, this cor- 
responds to heating that inflates the intracluster medium, 
driving the accretion shock to larger radii. 

Comparing the structure of shifted models to truncated 
models reveals interesting similarities. Figure 4 depicts 
the dimensionless density profiles for K c = 0.1, 0.3, and 
1.0 in clusters with c — 8. The run of density with radius 
is virtually identical for the shifted and truncated models 
when K c = 0.1 and is still quite similar when K c = 0.3. 
Thus, for most of the interesting range of entropy thresh- 
olds, there is little observable difference between shifted 
models and truncated models. However, as K c rises to 
1.0, the density profile of the shifted model becomes no- 
ticeably flatter than that of the truncated model. 




Fig. 4. — Density profiles for modified entropy models with c = 8 
and K c = 0.1, 0.3, and 1.0. Solid lines show /3(f) for truncated 
models (Kt), dotted lines represent shifted models (Kg), and the 
long- and short-dashed lines represent two types of radiative-loss 
model (J\r and Kqz). The density structure of truncated and 
shifted models is almost indistinguishable for small values of K c . 
The radiative-loss models have the elevated central densities char- 
acteristic of cooling- flow clusters. 

Similarity between the truncated and shifted models at 
the lower entropy thresholds arises because the unmodified 
entropy Kq depends approximately linearly on f g for small 
K. Thus, shifting the distribution by K c and truncating it 
at K c both produce similar modified entropy distributions 
and therefore similar density, temperature, and pressure 
distributions in the inner regions of the cluster. That is 
why entropy modification by heating and cooling can have 
similar effects on the luminosity and temperature of clus- 
ters as long as they establish similar values of the minimum 
entropy. 



2.3.3. Radiative Losses 

Observations show that many clusters contain gas with 
a cooling time less than a Hubble time (Fabian 1994), mo- 
tivating us to consider slightly more complicated modi- 
fications to the entropy distribution. Radiative cooling 
reduces the the specific entropy of a gas parcel according 
to 



ds_ _ d\nK 3 / 2 
~dl ~ dt 



lim p n p n e A(T) 

7r 



(16) 



This entropy equation can be rewritten to eliminate the 
dependence on gas density: 



dK 3 / 2 
dt 



TlpTl e 



{nm p )- 1 '' 2 T 1 l' 2 k{T) . (17) 



In general, the rate at which entropy declines will de- 
pend on the prevailing hydrodynamic conditions, stan- 
dard examples being isobaric cooling or isochoric cooling. 
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However, for the purposes of this toy model, let us as- 
sume that the right-hand side of this equation remains 
constant. We make this rather arbitrary assumption for 
two physical reasons. First, within the temperature range 
0.1 kcV < T < 2kcV, the quantity T 1 / 2 A(T) is approx- 
imately constant. Second, as gas settles within a nearly 
isothermal potential, work done by compression will at- 
tempt to maintain the temperature of that gas near the 
characteristic temperature of the halo. Thus, we are mo- 
tivated to use the following modification to represent ra- 
diative losses: 

K R (f g ) = {[k (f g + /c )] 3 / 2 - . (is) 

This modification corresponds to reducing the value of 
K 3 / 2 by the constant amount K^ 2 throughout the en- 
tropy distribution and discarding all the gas with K < 0. 
The short-dashed line in Figure 2 shows an example of the 
resulting distribution for c = 8 and K c = 0.3, and a similar 
line in Figure 3 shows how the maximum radius depends 
on K c . 

It will also be interesting to investigate a generalization 
of this radiative-loss model: 

k Ga (f a ) = {[Ko(f g + fcT - Kf ,a , (19) 

in which the tunable parameter a enables us to adjust the 
amount of gas below the cooling threshold. In the limit 
a — > co, this generalized distribution becomes identical to 
the truncated distribution. In what follows, we will fo- 
cus on the intermediate case of a = 3, represented by the 
symbol ATq3 and long-dashed lines in the relevant figures. 

Figure 4 compares the density distributions that re- 
sult from the radiative-loss models with the truncated and 
shifted models. For K c = 0.1, all the density profiles are 
very similar at f > 0.1, because the entropy profiles are 
nearly identical there. However, within that radius, en- 
tropy in the radiative-loss models drops below K Cl allow- 
ing for considerably greater compression and higher central 
density. This trend becomes more pronounced at higher 
k c because the proportion of the intracluster medium be- 
low the cooling threshold is even larger. 

2.4. Entropy Modification and Surface Brightness 

Modified-entropy models generated through shifts and 
truncations of the unmodified entropy distribution pro- 
duce surface-brightness profiles that are strikingly simi- 
lar to the /3-model profiles of many real clusters, at least 
within about 30% of the virial radius. Figure 5 shows the 
dimcnsionlcss emission measure S as a function of pro- 
jected radius f for three different truncated models, each 
with c = 8. As the threshold entropy progresses from 
k c = 0.1, characteristic of clusters, through k c — 0.3 
to k c — 1.0, characteristic of groups, the S(f±) profile 
substantially flattens. As long as A(T) does not change 
dramatically with radius, these emission-measure profiles 
will be nearly identical to the surface-brightness profiles. 

For comparision, Figure 5 also shows some /3-model pro- 
files selected to match the modified-entropy models. These 
/3-model profiles are defined by S oc [1 + (r/r c ) 2 ]~ 3/3+1 / 2 , 
where r c is the core radius and (3 determines the asymp- 
totic slope (Cavaliere & Fusco-Femiano 1978). In each 
case, the /3-model closely tracks the modified-entropy 
model within 0.3r 2 oo- Beyond that radius, the emission 
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Fig. 5. — Dimensionless emission measure S as a function of 
projected radius f±. The solid lines show S(f) for three modified 
entropy models with c = 8, truncated at K c = 0.1, 0.3, and 1.0. 
The dotted and dashed lines show emission measure profiles de- 
rived from /3-modcls with S oc [1 + (r /r c ) 2 ]- 3 f 3+0 - 5 . Notice that 
the /3-models are very good approximations to the modified-entropy 
models within 0.3r200- 

measures of the /3-models exceed those of the modified- 
entropy models, but owing to the low surface brightness at 
large radii, the /3-model parameters of clusters and groups 
have generally been measured within 0.3 — 0.5r2oo- Recent 
measurements by Vikhlinin, Forman, & Jones (1999) sug- 
gest that cluster surface brightness does indeed fall below 
the best fitting /3-model near the virial radius. 

In order to determine how the best-fitting (3 and r c 
depend on c and K Cl we performed surface-brightness 
weighted fits over a grid of modified entropy models. Be- 
cause the best-fit parameters depend somewhat on the por- 
tion of the profile being fitted, we elected to do these fits 
over two different radial intervals: 0.03 < f±^ < 0.3 and 
0.05 < fj_ < 1.0. Results for truncated (i^r) models 
are shown in Figure 6, and those for shifted (k$) mod- 
els are shown in Figure 7. Because the profiles in the 
shifted models tend to be shallower, the corresponding (3 
values are somewhat smaller, particularly for low values of 
k c . However, in both cases, the best-fitting (3 values span 
the observed range for clusters and groups of galaxies (see 
§ 3.2.1). 



2.5. Entropy Modification and Cluster Temperature 

One consequence of raising the core entropy of a clus- 
ter is a boost in the luminosity-weighted temperature of 
the intracluster gas. Figure 8 shows that increasing k c 
raises the temperature of gas at all radii in both the 
shifted and truncated models. Amplification of the cen- 
tral temperature can be quite extreme, exceeding 3T200 
when k c rises to 1.0. However, the resulting enhancement 
of the luminosity-weighted temperature is not nearly as 
dramatic. The value of T lum turns out to depend more 
critically on the concentration parameter c than on the 
entropy threshold, as demonstrated in Figure 9. 
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Fig. 6. — Best-fitting /3-modcl parameters for truncated modified- 
cntropy models. Each line gives the best-fitting (/?, r c ) pair for a 
particular value of the concentration c. The points on each line show 
best fits for five different values of the entropy threshold. The pro- 
gression K c = 1.0, 0.5, 0.3, 0.2, 0.1 runs from the upper left to lower 
right in each case. Because the best fits depend somewhat on the 
fitted interval, we depict best fits over the interval 0.03 < r j_ < 0.3 
with solid triangles and solid lines and best fits over the interval 
0.05 < fx < 1.0 with empty squares and dotted lines. 
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Fig. 8. — Dimensionless temperature profiles. Solid lines show 
how temperature (T) depends on radius [f ) in models with trun- 
cated entropy distributions (Kt)- Dotted lines show how tempera- 
ture depends on radius in models with shifted entropy distributions 
(Kg). Raising the entropy threshold K c increases the temperature 
at all radii for both kinds of entropy modification. This increase is 
slightly larger for shifted models because the extra gas at large radii 
leads to slightly higher central pressures. 
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Fig. 7. — Best-fitting /3-model parameters for shifted modified- 
cntropy models. Each line gives the best-fitting (/3, r c ) pair for a 
particular value of the concentration c. The points on each line show 
best fits for five different values of the entropy threshold. The pro- 
gression K c = 1.0, 0.5, 0.3, 0.2, 0.1 runs from the upper left to lower 
right in each case. Because the best fits depend somewhat on the 
fitted interval, we depict best fits over the interval 0.03 < f± < 0.3 
with solid triangles and solid lines and best fits over the interval 
0.05 < f± < 1.0 with empty squares and dotted lines. 



The reason that raising the core entropy of gas in an 
NFW potential has only a modest effect on its luminosity- 
weighted temperature involves both the gas density profile 



and the dark-matter density profile. Extra entropy flat- 
tens the density profile of intracluster gas, diminishing the 
relative luminosity of gas in the core. The characteristic 
radius of gas contributing the bulk of the cluster's lumi- 
nosity therefore moves outward. Because the shape of an 
NFW halo is steeper than isothermal at f > 1 /c, the char- 
acteristic temperature of gas within such a halo declines 
outside that radius (see Figure 8). Thus, the rise in gas 
temperature resulting from an increase in K c is largely 
offset by an increased contribution to the total luminos- 
ity from lower-temperature gas at larger radii. This ef- 
fect is most pronounced in the shifted models, which have 
substantial amounts of relatively cool gas beyond f = 1. 
In those models, Ti um actually declines as K c approaches 
unity. 

The large central temperatures in the shifted and trun- 
cated cases also seem somewhat unrealistic in light of ob- 
servations showing either level or radially increasing tem- 
perature gradients within cluster cores (e.g., Arnaud et al. 
2001a, b; Allen et al. 2001). In order to compare these tem- 
peratures more directly with observations, we must look 
at T± , the luminosity- weighted temperature along various 
lines of sight through a cluster. Figure 10 illustrates how 
T± depends on projected radius f± in clusters with c = 8, 
K c = 0.3, and four different kinds of entropy modification. 
Projected temperatures at f± < 0.1 in both the shifted 
and truncated models are still quite high (Tj_ w 2T2oo)- 
However, the core temperatures in the radiative-loss mod- 
els look more like those in real clusters (see § 3.2.2), rising 
with radius in the innermost regions and peaking between 
I.2T200 and I.6T200 at ~ 0.05r 2 oo in these particular mod- 
els. 
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Fig. 9. — Dimcnsionless luminosity- weighted temperatures (Ti um ) 
as a function of entropy threshold (K c ) for halos of concentration 
c = 4, 6, 8, 10, and 12. Solid lines depict truncated models (K^); 
dotted lines depict shifted models (Kg)- The dependence of Ti um on 
K c is rather modest over the interesting range 0.1 < K c < 1.0, es- 
pecially for larger values of c. This effect arises because a larger pro- 
portion of the emission from high-X c clusters comes from f > 1/c, 
where temperatures are lower than in the core. This outward shift 
of surface brightness largely offsets the rise in temperature at all 
radii evident in Figure 8. 
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Fig. 11. — Dimcnsionless luminosity-weighted temperatures for 
the standard radiative-loss model (Kb,)- Comparision of these 
Ti um (K c ) relations with those in Figure 9 shows that emission from 
the core gas lowers the luminosity- weighted temperature by ~30%, 
relative to the shifted and truncated models. 
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Fig. 10. — Dimcnsionless luminosity-weighted linc-of-sight tem- 
perature T± as a function of projected radius f±. The solid and 
dotted lines representing truncated and shifted models, respectively, 
both begin at rather high temperature (fa 2T200) at small radii and 
proceed to decline monotonically with radius. The long- and short- 
dashed lines representing radiative-loss models (see § 2.3.3) begin at 
low temperature at f = and rise to a maximum at 0.05r200, in 
better accord with observed clusters (see § 3.2.2). 
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Fig. 12. — Dimcnsionless luminosity-weighted temperatures for 
generalized radiative-loss model i\G3 with a = 3. Solid lines 
show the T\ UIn (K c ) relations for truncated models, and long-dashed 
lines show the corresponding relations for the radiative-loss models. 
Luminosity-weighted temperature is ~10% smaller in these gener- 
alized radiative-loss models. The difference is less dramatic than in 
the a = 3/2 case because the amount of gas below K c is smaller, 
resulting in somewhat higher core temperature and somewhat lower 
core density and luminosity. 
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Cool gas in the cores of clusters with entropy modi- 
fied through radiative losses can significantly affect the 
luminosity- weighted temperature if the fraction of gas be- 
low K c is relatively large. Figure 11 shows that T| U m 
for a given concentration c and entropy threshold K c is 
^30% lower in the standard radiative-loss model (-Kr) 
than in the corresponding truncated model. This effect 
is also present but not nearly so strong in the generalized 
radiative-loss model with a = 3, where T lum is <~10% lower 
than in the truncated models (see Figure 12). 



2.6. Entropy Modification and Cluster Luminosity 

All of our schemes for entropy modification substantially 
lower the luminosities of clusters. Figure 13 shows how di- 
mensionless luminosity L depends on the entropy thresh- 
old K c for clusters with c = 8. In all cases, luminosity 
declines by over an order of magnitude as K c rises from 
0.1 to 1. This decline is shallowest for the shifted model 
because the additional gas at large radii contributes sig- 
nificantly to the total luminosity when K c ~ 1. Note also 
that L is over twice as large in the standard radiative-loss 
case (Kr) because a large proportion of the luminosity 
in such models comes from the dense core. Luminosity 
enhancement in the a = 3 radiative-loss model is more 
moderate — only a few tens of percent greater than in the 
truncated case. 

Analyzing how L depends on K c leads to some impor- 
tant insights into how the cooling threshold determines the 
luminosity-temperature relation of clusters and groups. 
Figure 14 shows that L oc Kc 3 ^ 2 is a good approximation 
to the L-K c relation for truncated models in the range 
0.1 < K c < 1.0. This scaling reflects the asymptotic slope 
of the NFW density profile. Notice that the outer por- 
tions of the density profiles in Figure 1 track the underlying 
NFW profile of the dark matter, and the inner portions are 
nearly isentropic with K w K c . One can therefore approxi- 
mate the density profile of a truncated model by joining an 
isentropic density profile to an NFW profile at the radius 
?k at which the density of the unmodified NFW profile 
is pk ~ K c 3 ^ 2 . The analytical isentropic solution from 
§ 2.1.1 tells us that p ~ [f _1 ln(l + cf)] 3 / 2 within fx, so 
the luminosity from the isentropic gas scales oc p\r\. The 
luminosity from the NFW portion of the profile scales in 
the same way for the radii of interest. Thus, for larger val- 
ues of r , where the density of the NFW profile approaches 
p oc f~ 3 , we have L oc pK oc K c 3 ^ 2 - 

This particular scaling of L with K c has some very inter- 
esting consequences. Section 2.1.3 showed that the phys- 
ical luminosity of a cluster scales as L oc r\ w p\ r tvL . In 
order to substitute physical quantities for L, we need to 
recognize that K^ 2 oc K^ 2 T^/ 2 ' p cr and that K^ 2 oc 
T 1/,2 AtH, where tu is the age of the universe. We there- 
fore obtain 




Fig. 13. — Dependence of dimcnsionlcss luminosity L on the en- 
tropy threshold K c for clusters with c = 8. The decline of L with 
increasing K c is strong for all types of entropy modification. The ef- 
fect is stronger in the truncated (Kt) and radiative-loss (Xpj,, Kc:i) 
models than in the shifted model (Ks) because the latter has more 
gas at large radii, which contributes significantly to L when the 
entropy threshold is high. 
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Fig. 14. — Dimcnsionlcss luminosity L as a function of entropy 
threshold K c in truncated (-K"t) models for clusters with concen- 
tration c = 4, 6, 8, 10, and 12. The cluster luminosities in these 
models sharply decline as K c increases. For 0.1 < K c < 1.0 this 



decline approximately scales as L oc K, 



-3/2 



(20) 



where H is the Hubble constant at time in- 
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The following three features of this expression are par- 
ticularly noteworthy: 

• The power-law slope it predicts for the X-ray 
L oc T b relation is very close to the observed 
slope of b w 2.6 - 2.9 (e.g., Markevitch 1998; 
Arnaud & Evrard 1999). Because the halos of 
low-temperature groups have higher concentrations 
than high-temperature clusters (see § 3.1), we 
expect Ti um to decline by a factor ~ 1.3 as T\ um 
rises by an order of magnitude. This dependence 
on concentration steepens the power-law slope 
implied by equation (20) to b w 2.8. 

• The L-T relation does not depend on the form of 
the cooling function, which scales as A oc T 1 / 2 in 
the free-free cooling regime at T > 2 keV and as 
A oc T -1 / 2 in the line-cooling regime at T < 2kcV. 
If the entropy threshold in real clusters did not 
depend on A, then we would expect a distinct 
steepening of the L-T relation as T rises through 
<~ 2keV, which is not observed. Equation (20) 

is indepedent of A because the dimensionless 
gas mass within tk scales oc Pkt\, and thus 
remains approximately constant in an NFW halo 
for fi > 1 /c. The cooling time of that gas is <~ tn, 

so its physical luminosity scales ~ M 2 ooT2oot : a 1 ■ 
This scaling law also applies to the entire cluster, 
because the luminosity of gas outside tk scales in 
lockstep with gas inside tk ■ Applying the relation 
M 2 oo oc T^qqH -1 thus leads back to equation (20), 
except for the T\ um factor. 

• Finally, the L-T relation implied by equation (20) 
changes very little with time, again in agreement 
with observations (Mushotzky & Scharf 1997; 
Donahue et al. 1999; Delia Ceca et al. 2000; 
Borgani et al. 2001b). That happens because the 
mass of a cluster of a given temperature declines 
oc H^ 1 with increasing redshift. According to the 
L oc -A^ooT^ooifj 1 scann S) the resulting decline 

in the intracluster gas mass compensates almost 
precisely for the drop in the cooling threshold 
owing to the shorter cooling time tu- 

This scaling breaks down at the low-temperature end, 
where K c > 1, because the entire intracluster (or in- 
tragroup) medium in the truncated and shifted cases 
becomes effectively isentropic. Then, the gas den- 
sity determined by the cooling threshold is oc T/Atu, 
and the luminosity from within r 2 oo scales as L oc 

r i™ A ~ 2 ^ /2 (#tH)~ 2 # -1 - ThuS > WC Sh ° uld CX P CCt th ° 

L-T relation for low-temperature halos to be steeper than 
for clusters, roughly L oc T 4 - 5 for A ~ T -1 / 2 , and we 
should expect the luminosity of objects at a given temper- 
ature to increase with time. 

In order to make further progress, we will need to un- 
derstand how halo concentration depends on halo mass. 

3. MODELING REAL CLUSTERS 

The previous section explored the properties of dimen- 
sionless modified-entropy models for clusters depending 



on two parameters, the halo concentration parameter c 
and an entropy-threshold parameter K c . As in Voit & 
Bryan (2001), we argued that the crucial entropy thresh- 
old should be the entropy level at which intracluster gas 
would cool in a Hubble time. Now we will use a set of re- 
lations between halo concentration and halo mass to con- 
struct modified-entropy models for real clusters. In these 
models, the halo mass M 2 qq determines the concentration 
parameter c. The mass and halo concentration then set the 
temperature of the cluster, and this temperature sets the 
level of the entropy threshold K c . Cluster properties are 
therefore determined primarily by the overall cosmological 
model, because it fixes the global baryon fraction, the rela- 
tion between halo mass and concentration, and the Hubble 
time that governs K C {T). Some of a cluster's properties, 
such as the temperature of the innermost gas, also depend 
on the chosen scheme for entropy modification. The prop- 
erties of the resulting modified-entropy models turn out 
to be remarkably similar to those of observed clusters and 
groups. 

The section begins by defining the M200-C relations we 
will use and comparing them to some of the scarce data. 
Then, we explore the structure of the intracluster medium 
in our modified-entropy cluster models, showing that the 
observed relationship between (3 and Ti um may be more 
than just a selection effect and demonstrating how en- 
tropy modification can lead to temperature profiles similar 
to those observed in both the inner and outer regions of 
clusters. We go on to compare the mass-temperature re- 
lation from our modified-entropy models to observations, 
showing that the models reproduce the observational re- 
sults of Nevalainen et al. (2000) and Finogucnov et al. 
(2001). We also compare the Lx-T\ um relation from our 
models to cluster data, again showing excellent agreement. 
Much of the dispersion of this relation, in both the models 
and the data, stems from differences in the amount of gas 
below the cooling threshold, and we conclude with some 
speculations about what governs these differences. 

3.1. Concentration Parameter 

The concentration of a dark matter halo is determined 
by its formation history. Early formation generally leads to 
a denser core and thus a higher concentration. Low-mass 
halos, which tend to collapse earlier in time, are therefore 
expected to be somewhat more concentrated than high- 
mass halos. Because a halo's collapse and merger history 
determines its concentration, the relation between concen- 
tration and halo mass depends on the underlying cosmo- 
logical model, specified by the matter density fijui the 
dark-energy density Q\, the normalization as of the per- 
turbation spectrum, and the shape of that spectrum, often 
parameterized by T (Bardcen ct al. 1986; Sugiyama 1995). 
In this paper, we will restrict our attention to currently- 
favored ACDM models, in which the resulting clusters turn 
out to be thoroughly consistent with observations. 

Several versions of the concentration parameter ca can 
be found in the literature, each defined with respect to 
a radius ta within which the mean halo density is Ap cr . 
Here we will define c with respect to r 2 oo, in order to be 
consistent with the models developed in the previous sec- 
tion; thus, c = C2oo- Figure 15 shows various relations 
between halo mass and C200 that can be found in the lit- 
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Fig. 15. — Concentration parameter (C200) as a function of 
halo mass (M20o)- Solid triangles give concentrations derived 
from Chandra observations of real clusters by Allen ct al. (2001). 
Empty squares give the concentrations of clusters simulated by Eke, 
Navarro, & Frenk (1998) in a ACDM cosmology with <j% = 1.1. Solid 
lines show C200 values predicted by models 50.9, Si. 2, and 51.6 from 
Eke et al. (2001), for a ACDM cosmology in which ct 8 = 0.9, 1.2, 
and 1.6, respectively, for a power spectrum with shape parameter 
T = 0.2. Dashed lines show C200 given by models To.i and To. 5 
from Eke ct al. (2001), for which V = 0.1 and 0.5, respectively, and 
erg = 0.9. The dotted line labeled TN shows the C200 values from 
the fit of Tozzi & Norman (2001) to the NFW prescription for a 
ACDM cosmology with erg = 1.1. 

erature. The Chandra observations of Allen et al. (2001) 
suggest that concentrations corresponding to C200 ~ 4 — 6 
are appropriate for massive (~ 10 15 M Q ) clusters. Concen- 
trations of groups are less certain because they are more 
model-dependent, and the dispersion in C200 is expected 
to increase at lower halo masses (Afshordi & Ccn 2002). 

In this section, we will adopt the halo mass- 
concentration relations predicted by models So. 9 and Si. 2 
from Eke, Navarro, & Stcinmetz (2001), which are con- 
sistent with the Chandra data on massive clusters. Both 
models are derived from A-body simulations of ACDM 
cosmologies with Qm = 0.3, ft a = 0.7, and h — 0.65. 
The power spectrum in each model has a shape parame- 
ter T = 0.2, but the normalization of the power spectrum 
differs: erg = 0.9 in model So. 9 and as = 1.2 in model Si. 2- 
The higher normalization in model Si. 2 leads to halos that 
are about 30% more concentrated. Figure 16 shows the de- 
pendence of K c on halo mass that follows from model Si. 2; 
the entropy thresholds for model So. 9 are very similar. 



3.2. Cluster Structure 

The cluster models we analyze in this section are com- 
pletely determined by the assumed relationship between 
halo mass and concentration (either So. 9 or Si. 2) and the 
chosen scheme for entropy modification (At, As, ATr, 
or Kqs). For each halo mass, we compute the bolomet- 
ric X-ray luminosity Lx, luminosity- weighted temperature 
Ti um , X-ray surface-brightness profile Sx(r±), and pro- 
jected temperature profile T±(r±). Unlike in § 2.1.3, each 
integral includes a cooling function A(T) drawn from the 
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Fig. 16. — Dimcnsionless entropy threshold (K c ) determined 
by cooling as a function of halo mass M200. The lines show the 
entropy thresholds for truncated models (solid line), shifted mod- 
els (dotted line), and radiative-loss models with a = 3/2 (short- 
dashed line), and a = 3 (long-dashed line), assuming the halo mass- 
concentration relation from model 51.2. The corresponding thresh- 
olds for model So. 9 are very similar. Above 10 14 /i~ 1 Mq, differ- 
ent entropy-modification schemes produce slightly different dimen- 
sionlcss thresholds because the temperature Ti um associated with a 
given halo mass is slightly different. Equations (12) and (13) give 
the entropy threshold K c and dimcnsionless entropy threshold K c , 
respectively, as functions of temperature. 

models of Raymond & Smith (1977) for gas with 30% 
solar metallicity. We will first focus on the spatially re- 
solved properties of these model clusters, showing how 
their surface-brightness profiles and temperature gradients 
compare with those of observed clusters. 

3.2.1. Surface- Brightness Profiles 

Modificd-cntropy models of groups tend to have flatter 
surface-brightness profiles than hot clusters, in accordance 
with observations (e.g., Ponman et al. 1999; Horner et 
al. 1999). Because the surface-brightness profiles of both 
real clusters and our modified-entropy models are well de- 
scribed by /3-models with Sx oc [1 + (r/r c ) 2 ]~ 3/3+1 / 2 inside 
of 0.3r2oo (see § 2.4), we compare our models with the 
data by comparing the best-fitting /3-model parameters. 
Figure 17 shows how [3 and r c depend on Ti um in trun- 
cated modified-entropy models with relation Si. 2 and in 
the cluster data of Finoguenov et al. (2001). The models 
reproduce the observed range of both parameters and the 
observed rise in each parameter with increasing tempera- 
ture. However, the observed dispersion in each parameter 
at a given temperature is quite large. Wu & Xue (2002) 
find the same kind of temperature dependence in their 
cluster models, based on the work of Bryan (2000), which 
are very similar to the truncated models developed in this 
paper. 

Some of the dispersion in [5 and r c at a given tempera- 
ture may stem from selection effects. Because the surface- 
brightness profiles of modified-entropy clusters continually 
steepen with radius, unlike those of /3-models, the best- 
fitting f3 and r c depend somewhat on the range of radii 
being fit. Figure 17 shows that the best-fitting values of [3 
and r c both decline as the surface-brightness threshold of 
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Fig. 17. — Best-fitting /3-model parameters as a function of 
luminosity- weighted temperature Ti um . In a /3- model, X-ray sur- 
face brightness declines with radius as Sx oc [1 + (r/r c ) 2 ] _3 ' 3+1 ' /2 , 
where r c is the core radius and /3 determines the asymptotic slope. 
Triangles show measured values of /3 (upper panel) and r c (lower 
panel) compiled by Finogucnov ct al. (2001). Solid lines show 
the values of these parameters in the best fits to truncated mod- 
ified entropy models (Kt) over the radial range 0.03 < f < 0.3, 
where f = r/r2oo • Dotted lines show the best-fitting values over the 
range 0.05 < f < 1.0. The other lines show best fits above surface- 
brightness thresholds of 2 X 10~ 15 ergcm~ 2 s -1 arcmin" 2 (short- 
dashed lines), 2 X 10 -14 erg cm -2 s _1 arcmin -2 (long-dashed lines), 
and 2 X 10 -13 erg cm -2 s -1 arcmin -2 (dot-dashed lines). 
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Fig. 18. — Best-fitting /3-model parameters to various modified 
entropy models as a function of luminosity-weighted temperature 
Tlum- As in Figure 17, triangles show measured values of f3 (upper 
panel) and r c (lower panel) compiled by Finogucnov et al. (2001). 
Solid lines show the values of these parameters in the best fits 
to truncated modificd-entropy models (Kt), dotted lines show the 
best-fitting values for shifted models (Kg), short-dashed lines show 
best fits for radiative- loss model K n, and long-dashed lines show 
best fits for radiative-loss model Kq^. All fits were performed over 
the radial range 0.03 < f < 0.3 for models based on concentration- 
parameter relation Si. 2- 



the fit rises from 2 x 1CP 15 erg cm" 2 s" 1 arcmin -2 through 
2 x 10~ 14 ergcm~ 2 s" 1 arcmin -2 to 2 x 10~ 13 ergcm~ 2 s" 1 
arcmin -2 . This effect is more severe for groups than for 
clusters because their overall surface brightness is smaller. 
The figure also shows that fits over the radial range 0.05 < 
r/r2oo < 1.0 yield slightly larger (3 and r c than fits over 
0.03 < r/r 200 < 0.3. Komatsu & Seljak (2001) noted 
a similar surface-brightness selection effect in their poly- 
tropic models for clusters. However, the rise of (3 with T\ um 
in their models is entirely due to this surface-brightness 
bias. In our models, it is largely a consequence of the 
greater impact of entropy modification on groups. Puz- 
zlingly, this selection effect becomes stronger at low sur- 
face brightness levels in the models of Wu & Xue (2002), 
perhaps because their boundary condition leads to a less 
precipitous decline in Sx near the virial radius. 



The observed dispersion in j3 and r c may also depend 
in part on the amount of gas below the cooling threshold. 
Figure 18 illustrates how (3 and r c depend on the scheme 
for entropy modification. Shifted models have core radii 
identical to those of the corresponding truncated models 
and slightly smaller values of (3. However, allowing for 
radiative losses reduces the best-fitting core radius by a 
factor of 3 — 5 and flattens the (3-Ti um relation, yielding 
w 0.5 - 0.6 for fits over the interval 0.03 < r/r 20 o < 0.3. 
Other possible sources of dispersion include ellipticity and 
cluster substructure, which cannot be investigated in the 
context of our spherically symmetric models. 
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Fig. 19. — Luminosity-weighted projected temperature T± as a 
function of projected radius r±. Radii are given in units of the radius 
f*2500 within which the mean matter density is 2500 times the critical 
density. Temperatures are given in units of the mean mass-weighted 
temperature T2500 within T2500- The hatched area shows the best fit 
of Allen et al. (2001) to their Chandra data on six hot clusters with 
strong cooling flows, along with the uncertainty in that fit. Other 
lines show Tj_ from modified-entropy models with parameters char- 
acteristic of hot clusters: c = 4, K c = 0.1, and T200 = lOkeV. The 
short-dashed line depicts a standard radiative-loss model (Kr), the 
long-dashed line depicts generalized radiative-loss model Kg3 with 
o = 3, the dot-dashed line shows a generalized radiative-loss model 
with a = 6, and the solid line shows a truncated model (Kt), in 
which a is effectively infinite. 



3.2.2. Temperature Gradient at Small Radii 

The temperature gradient within ~ 0.3 r 2 oo also appears 
to depend on the amount of gas below the cooling thresh- 
old. Figure 19, in which the model parameters are ap- 
propriate for hot clusters, shows that T±(r±) gradually 
declines with radius when the entropy distribution is sim- 
ply truncated at the entropy threshold. Our radiative- 
loss models, in contrast, have temperatures that rise out- 
ward from r = 0, then flatten, becoming virtually isother- 
mal out to 30% of the virial radius. As the parameter 
a drops from the very large values characteristic of trun- 
cated models to the canonical radiative-loss value of 3/2, 
the region with a positive temperature gradient becomes 
progressively larger. The run of temperature with radius 
in our model with a — 3/2 is very similar to the "univer- 
sal" temperature profile observed by Allen et al. (2001) in 
clusters that appear to have strong cooling flows. In or- 
der to facilitate comparison of our models with that data, 
Figure 19 shows radii in units of r25uo 5 the radius within 
which the mean matter density is 2500p cr , and tempera- 
ture in units of T2500, the mean mass- weighted tempera- 
ture of intracluster gas within r25oo- For clusters of this 
concentration, r25oo ~ 0.3r2oo- Note that setting a = 6 
produces a cluster model with a temperature gradient re- 
maining within 10% of isothermal from small radii to be- 
yond 0.3r 20 o- 
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Fig. 20. — Effective polytropic index 7 g = din P/d\n p as a 
function of radius in modified-entropy models. Solid lines depict 
truncated models (A"t), an d dashed lines depict radiative-loss mod- 
els (Kc:j) with a = 3. Results are given for models of clusters with 
three different luminosity-weighted temperatures (Ti um = 0.5, 2.0, 
and 10.0 keV, top-to-bottom), using concentration-parameter rela- 
tion Si. 2- The curves show that the truncated models are iscntropic 
(7eff = 5/3) at small radii and nearly isothermal (7 e ff ~ 1-1 — 1-2) 
at large radii. Values of 7 e ff in the radiative-loss models are simi- 
lar at large radii but drop below unity in the radiative-loss models 
because of the low central temperatures. 



3.2.3. Temperature Gradient at Large Radii 

Projected cluster temperatures at large radii gradually 
decline with radius in all of our modified-entropy mod- 
els. ASCA observations of many clusters show evidence 
for such a decline, which is often expressed as an effective 
polytropic index 7 e ff = d\nP/d\np (e.g., Markevitch et 
al. 1998; Markevitch et al. 1999; Ettori & Fabian 1999). 
These observations indicate that 7 e ff ~ 1-2 outside the 
cores of clusters. However, recent XMM-Newton observa- 
tions have shown that at least some clusters are nearly 
isothermal within 0.4r2oo, except for the low-temperature 
core (Arnaud et al. 2001a, b). The results of § 3.2.2 suggest 
that radiative losses may be responsible for the isothermal 
profiles found in these clusters. 

Figure 20 shows how 7 e g varies with radius in truncated 
modified-entropy models and in radiative- loss model Kq^. 
The cores of the truncated models are nearly isentropic, 
with 7 ff w 5/3 at r < 0.01. This isentropic region extends 
to larger radii in cooler halos because their entropy thresh- 
olds amount to a larger fraction of the halo's mean entropy. 
Outside the isentropic region, we find 7 fr w 1.1 — 1.2, in 
agreement with observations. The origin of this value of 
7 ff can be traced back to our assumption that the unmod- 
ified gas and dark-matter density profiles are described by 
identical NFW laws. Komatsu & Seljak (2001) have ex- 
ploited that same assumption to construct polytropic mod- 
els for the intracluster medium, finding that 7 e fi ~ 1.1 — 1.2 
produces the best agreement between these density profiles 
at large radii. 

The behavior of j c g in the radiative-loss models is sim- 
ilar to the truncated models at large radii but quite dif- 
ferent at small radii. Figure 20 shows that the effective 
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polytropic index steadily rises from 7 e ff < 1 inside the core 
through 7 e ff w 1 near the core radius to 7 e fi ~ 1.1 — 1.2 at 
larger radii. The point at which a cluster makes this tran- 
sition from an isothermal gradient to a mildly declining 
gradient depends somewhat on cluster temperature. Hot 
clusters with massive halos tend to have smaller concen- 
trations which allow for near-isothermality at larger radii. 

In order to compare our models with observations, we 
have evaluated -f c s at r/r 2 oo — 0.2, essentially twice the 
core radius. This radius is large enough to escape the 
severe mo del- dependence of 7 e fj within the cluster core 
but small enough that the signal-to-noise in cluster ob- 
servations still allows for accurate measurements. Fig- 
ure 21 compares those evaluations of 7 e g- with the data 
of Finoguenov ct al. (2001). Our models agree with the 
data, within the large observational uncertainties, and also 
share the observed (2.5<j) tendency for 7 c ff to decline with 
increasing temperature. Two different effects combine to 
produce higher values of 7 e ff at low temperature. Because 
the halos of groups tend to be more concentrated than 
those of clusters, the scaled temperature T/T 2 oo tends to 
be larger in the cores of groups, leading to higher 7 e fj. An 
additional enhancement of j c g appears in truncated mod- 
els, in which the entropy threshold has a more significant 
effect at 0.2r2oo on group scales (see Figure 20). 




T lum (keV) 

Fig. 21. — Effective polytropic index 7 e ff as a function of 
luminosity- weighted temperature Ti um . The triangles show 7 e g val- 
ues derived from observed clusters by Finoguenov et al. (2001). The 
lines show values of 7 e // at r/r200 = 0.2 in truncated modificd- 

entropy models (Kt, solid line) and radiative- loss models (Kq3, 
dashed line). 



3.3. Mass -Temperature Relation 

Many efforts to measure cosmological parameters such 
as VLm and as have relied heavily on the cluster mass- 
temperature relation to link cluster masses to observables 
like cluster temperature and luminosity (e.g., Henry 1997; 
Bahcall & Fan 1998; Eke et al. 1998; Donahue & Voit 
1999; Borgani et al. 1999; Henry 2000). Until very re- 
cently, observers had to rely on M-T relations that were 
calibrated with numerical simulations (e.g., Evrard et al. 
1996; Bryan & Norman 1998). These simulations provide 
the constant of proportionality Ti um in the scaling relation 
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Values of T\ um from simulations in the literature can differ 
by as much as 30%; examples include T lum « 0.9 from 
Evrard et al. (1996), T| um w 0.8 from Bryan & Nor- 
man (1998), f lum w 0.9 from Thomas et al. (2001), and 
2] um ~ 0.6 from the non-radiative simulations of Muan- 
wong et al. (2001). The vast majority of these sim- 
ulations ignore radiative cooling and supernova heating, 
both of which can elevate Tium- Mass-temperature rela- 
tions derived from spatially-resolved X-ray observations 
now suggest that the normalization constant T\ um derived 
from simulations is indeed too low (Horner et al. 1999; 
Nevalainen et al. 2000; Finoguenov et al. 2001). Further- 
more, cluster-formation simulations that include radiative 
cooling also yield a higher normalization constant; for ex- 
ample, Ti um 1.1 in the radiative simulations of Muan- 
wong et al. (2001). 

Our modificd-cntropy models support these findings. 
Figure 22 compares two M5oo-Ti um relations derived from 
our modificd-entropy models with the data from Finoge- 
unov et al. (2001) and Nevalainen et al. (2000). In these 
relations, we scale to M 50 o, the halo mass inside the radius 
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Fig. 22. — Relation between M500 and luminosity- weighted 
temperature (T lum ). Solid triangles show mass measurements by 
Finoguenov ctal (2001), using /3-modcl fitting and the assumption 
of hydrostatic equilibrium. Open squares show mass measurements 
by Nevalainen et al. (2000), also using /3-model fitting and hydro- 
static equilibrium. Both of these data sets were corrected for the 
presence of cooling flows, so wc compare them to our models without 
radiative losses. The solid and dotted lines show the M5oo-7i um re- 
lations derived from Kt models using the concentration parameter 
relations Si. 2 and So. 9, respectively (see Figure 15). The corre- 
sponding relations derived from shifted models (Kg) are virtually 
identical. The dashed line shows the M5oo-7i um relation derived by 
Evrard et al. (1996) from simulations that do not include radiative 
cooling or supernova feedback. 
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^500 within which the mean matter density is 500p cr , be- 
cause measurements of M 500 are thought to be more re- 
liable than measurements of M 2 oo- The data have been 
corrected for any cooling flow that might be present in the 
cluster core, so we compare them to our truncated models 
with concentration laws So. 9 and Si. 2- The correspond- 
ing M 5 oo-Tl um relations generated from shifted models are 
virtually identical. Both the truncated and shifted models 
agree well with the data, but those based on concentration 
relation Si. 2 appear to match the data slightly better. 

This agreement between the observed M 5 oo-7i um rela- 
tion and the modified-entropy models is quite compelling, 
considering that the models contain no adjustable param- 
eters. For comparision, Figure 22 also shows the relation 
derived from simulations by Evrard et al. (1996), which 
has a lower temperature normalization (equivalent to a 
higher mass normalization) and a shallower slope. The 
steeper slope of the M^qq-T relation from modified-entropy 
models arises primarily from the dependence of T\ um on 
halo concentration. The higher concentration of low-mass 
halos leads to a higher luminosity-weighted temperature 
for a given mass. Komatsu & Seljak (2001) find a similar 
effect in their polytropic cluster models. Entropy modi- 
fication produces some additional steeping in the relation 
and shifts its normalization to higher temperatures (see 
Figure 9). 

The greater dispersion and offset toward lower masses 
at Ti um ?y 1 seen in Figure 22 could be due to surface- 
brightness bias (Lloyd-Davies, Bower, & Ponman 2002). 
Because the gas temperature of intracluster gas in our 
models tends to decline with radius at large radii, the 
luminosity-weighted temperature depends on the maxi- 
mum radius over which that temperature is measured. 
This effect is negligible in clusters because the bulk of the 
emission is concentrated around the cluster's core, but it 
can be significant in groups. Figure 23 shows how the 
observed M 500 -Ti um relation would depend on the limit- 
ing surface-brightness for truncated models with Si. 2 con- 
centrations. Increasing the surface brightness threshold 
changes the derived relation in the sense that is observed. 
Figure 24 shows the relations for different entropy modi- 
fication schemes at a fixed surface-brightness threshold of 
5 x 10~ 15 ergcm~ 2 s _1 arcmin~ 2 . Truncated and shifted 
models are affected similarly, but surface-brightness bias 
has less effect on the radiative-loss model (Kg3) because 
its central temperature gradient is closer to isothermal. 
Some of the scatter to higher temperatures may also arise 
from the scatter in concentrations expected in low-mass 
halos (e.g., Afshordi & Cen 2002), but the analysis of 
Lloyd-Davies et al. (2002) suggests that halo concentra- 
tion alone cannot account for the large dispersion of the 
mass-temperature relation for groups. 

Despite the excellent agreement between our models and 
the observations, a word of caution is in order. Both the 
modified-entropy models and the X-ray mass derivations 
assume that clusters are in hydrostatic equilibrium. Sim- 
ulations, on the other hand, suggest the intracluster ve- 
locity field is not completely relaxed, with turbulent ve- 
locities ~ 10 — 20% of the sound speed (Packer & Sarazin 
2001). Furthermore, the luminosity-weighted temperature 
derived from theoretical modeling is not necessarily iden- 
tical to the spectral-fit temperature derived from obser- 



vations (Mathiesen & Evrard 2001). Nevertheless, our 
results strongly imply that entropy modification and the 
cosmology-dependent relationship between halo mass and 
concentration must be accounted for in simulations de- 
signed to normalize the cluster mass-temperature relation. 
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Fig. 23. — Relation between M500 and luminosity- weighted tem- 
perature (Ti um ) at different surface-brightness levels. Because tem- 
perature declines with radius at large distances from the center of 
our group models, the limiting surface brightness affects the ob- 
served mass-temperature relation. The solid line shows the M500- 
Ti um relation for model Kt,S\.2- The other lines show how that 
relation changes when temperature is measured within regions where 
the bolomctric surface brightness in ergem - 2 s _1 arcmin - 2 exceeds 
1 X lCT 15 (dotted), 2 X lO" 15 (short-dashed), and 5 X 1CT 15 (long- 
dashed). The data are the same as in Figure 22. Some of the ob- 
served dispersion below the solid line could stem from this surface- 
brightness bias. 




Fig. 24. — Relation between M500 and luminosity- 
weighted temperature (Ti um ) above a surface brightness of 5 X 
10~ 15 crgcm~ 2 s _1 arcmin -2 . The solid and dotted lines show the 
relation for truncated and shifted models, respectively. The dashed 
line shows the relation for radiative- loss model Kgz- The data are 
the same as in Figure 22. The presence of gas below the cool- 
ing threshold near the center of the radiative-loss model lowers the 
central temperature, thereby mitigating the effects of the surface- 
brightness threshold. 
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3.4. Luminosity- Temperature Relation 

Our previous work on modified-entropy models and the 
cooling threshold demonstrated how well these models ac- 
count for the observed slope and normalization of the clus- 
ter luminosity-temperature relation (Voit & Bryan 2001). 
That paper compared truncated and shifted modified- 
entropy models to cluster data that were corrected for the 
effects of cooling flows. Here, we revisit those models and 
show how the Lx~Ti um relation changes when gas below 
the cooling threshold is included. 

Figure 25 shows the Lx-T\ um relation for clusters with- 
out cooling flows. Of the three models shown, the trun- 
cated model with concentration relation Si. 2 best fits the 
data, but all three models faithfully describe the Lx-T\ unl 
relation equally well on cluster scales. Substituting con- 
centration relation So. 9 shifts the models to slightly lower 
temperatures, but the agreement with groups is almost as 
good. Using shifted instead of truncated entropy distribu- 
tions elevates the X-ray luminosity on group scales by a 
factor <~ 2. 

While there arc essentially no free parameters in these 
models, the cluster luminosity does depend on our as- 
sumed value of Qm ■ Raising £1 m would lower the global 
ratio of baryons to dark matter, thereby lowering the lumi- 
nosities of clusters in our models oc flj^ (see Eq. 9). Con- 
versely, lowering Qm would raise the cluster luminosities. 
The excellent agreement found here for Am = 0.33 should 
therefore reinforce confidence in similar measurements of 
the matter density from the baryon-to-dark-matter ra- 
tios of clusters (e.g., Evrard 1997). Our models imply 
that groups are less suitable for such measurements be- 
cause their baryon-to-dark-matter ratio within r2oo may 
be < 50% of the global value (see Figure 3). 

The tendency of the data at ~ 1 keV to scatter toward 
the low-Lx, high-Xi um side of the model relations may be 
due to the same effects that cause similar scatter at those 
temperatures in the observed M 50 q-Tx relation. Figure 26 
shows how surface-brightness bias affects the luminosity- 
temperature relations derived from our models. Raising 
the surface-brightness threshold within which luminosity 
and temperature are measured lowers Lx while raising 
Tium, producing the same trend seen in the data. Mea- 
suring both of these quantities only within regions where 
Sx > 5 x 10~ 15 ergcm~ 2 s _1 arcmin -2 produces relations 
that track the lower envelope of the group data. Unlike the 
cluster measurements, the data on groups have not been 
corrected for the effects of gas below the cooling threshold, 
which could counterbalance the surface-brightness bias by 
raising Lx and lowering Ti um (see Figure 27). We there- 
fore suspect that the dispersion in the L-T data of Hclsdon 
& Ponman (2000) arises from a combination of differing 
surface-brightness biases and differing amounts of gas be- 
low the cooling threshold. 

Reproducing the behavior of clusters and groups with 
significant amounts of gas below the cooling threshold re- 
quires an additional parameter. In our models, that pa- 
rameter is a, which characterizes the slope of the entropy 
distribution function below the cooling threshold. Fig- 
ure 28 shows the relationship between Ti um and £rosat, 
the ROSAT-band (0.1-2.4 keV) luminosity, in the data 
from Markevitch (1998) and in models with differing val- 
ues of a. Triangles show data corrected to remove the 
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Fig. 25. — Relation between bolometric X-ray luminosity Lx 
and luminosity- weighted temperature (Ti um ). Solid triangles show 
measurements of clusters with insignificant cooling flows compiled 
by Arnaud & Evrard (1999). Open squares show cooling-flow cor- 
rected measurements by Markevitch (1998). Solid circles show group 
data from Helsdon & Ponman (2000). Because the cluster data are 
not strongly influenced by cooling flows, we compare them to mod- 
els without radiative losses. The solid and dotted lines show the 
Lx-T\ um relations derived from truncated (Kt) and shifted (Ks) 
models, respectively, and the concentration parameter relation Si. 2 
(see Figure 15). The dot-dashed line shows the relation derived 
from model Kt using the concentration parameter relation So. 9. 
The models assume a standard ACDM cosmology with Qm = 0.33, 
SIa = 0.67, and Qf, = 0.02 h -2 , and a Hubble parameter of h = 0.65 
has been applied to both the models and the data. 
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Fig. 26. — Relation between Lx and luminosity-weighted tem- 
perature (T lum ) at different surface-brightness levels. Because tem- 
perature declines with radius at large distances from the center of 
our group models, the limiting surface brightness affects the ob- 
served mass-temperature relation. The solid line shows the Lx- 
Ti um relation for model Kt,S\.2- The other lines show how 
that relation changes when temperature and luminosity are mea- 
sured within regions where the bolometric surface brightness in 
ergcm~ 2 s~ 1 arcmirr " 2 exceeds lxl0~ 15 (dotted), 2xl0~ 15 (short- 
dashed), and 5 X 10~ 15 (long-dashed). The data are the same as 
in Figure 25. Some of the observed dispersion below the solid line 
could stem from this surface-brightness bias. 



18 



MODIFIED-ENTROPY MODELS 









1 








= 0.33 








- 


"A 


= 0.67 










n b 

h 


= 0.02 h- 2 

= 65 --^ ^ 


■J 


□ * 














; 


- 
E 

1 


/ 








: 
: 

; 


z 




if 'i*/ S •> S x 1 f 1 


-15 

K 7 




= 




H 






°1.2 












S 1.2 








i 


• ^1.2 
1 








i 




10 








T lum (keV) 









Fig. 27. — Relation between Lx and luminosity-weighted 
temperature (Ti um ) above a surface-brightness threshold of 5 X 
10~ 15 crgcm~ 2 s _1 arcmin~ 2 . The solid and dotted lines show the 
relation for truncated and shifted models, respectively. The short- 
dashed and long-dashed lines show the relations for radiative-loss 
models and Kgs, respectively. The data are the same as in 
Figure 25. The presence of gas below the cooling threshold near the 
centers of the radiative-loss models lowers the central temperature 
and raises the central luminosity. Much of the dispersion to the up- 
per left of the solid and dotted lines near ~ 1 keV could be due to 
cooling flows, which have not been corrected for in the group data. 
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Fig. 28.— Relation between ROSAT (0.1 - 2.4 keV) X-ray lumi- 
nosity Lrosat and luminosity- weighted temperature (Ti um ). Open 
squares show uncorrected measurements from Markevitch (1998), 
and solid triangles connected to those squares show results for those 
same clusters after correcting for a central cooling flow. (Error 
bars representing temperature uncertainty ~ 0.5 keV have been sup- 
pressed for legibility.) The solid line shows the I/ROSAT-Tlum rela- 
tion for truncated models, the long-dashed line depicts model Kg3 
(a = 3), and the short-dashed line depicts model Kr {a = 3/2). 
Dotted lines connect models for identical halo masses, showing how 
the iR,oSAT-Ti um relation for a given halo mass shifts to higher 
luminosity and lower temperature as a declines. 



effects of gas below the cooling threshold, and the squares 
connected to those triangles show the uncorrected data. 
While the corrected data cluster around the line repre- 
senting truncated models, many of the uncorrected data 
points lie closer to the line representing radiative-loss mod- 
els with a — 3. A few uncorrected data points stray even 
farther from the truncated models, toward models with 
a < 3. 

Notice that the lines connecting data points with the 
largest corrections are generally parallel to the lines con- 
necting modified-entropy models with identical halo mass. 
This agreement suggests that the parameter a adequately 
characterizes the offset in Lx-T\ um space owing to intr- 
acluster gas below the cooling threshold. An interesting 
test of these radiative-loss models would be to compare 
the value of a inferred from the luminosity-temperature 
offset to that implied by the inner temperature gradient. 

3.5. What governs a? 

The modified-entropy models we have constructed sug- 
gest that most of the observable properties of clusters de- 
pend on only two parameters. The halo mass M200 deter- 
mines a sort of cluster "main sequence" with well-defined 
-^500-Tlum and Lx-T\ um relations. Deviations from those 
relations and the severity of the inner temperature gradi- 
ent both depend on a second parameter, a, related to the 
amount of intracluster gas below the cooling threshold, 
but what physical processes determine the value of a? 

One is tempted to interpret a in terms of a mass cooling 
rate Mx', however, XMM-Newton observations show little 
evidence for the low-temperature line emission expected in 
the standard cooling- flow picture (Peterson et al. 2001). 
Thus, instead of viewing a as simply a measure of cooling- 
flow strength, we would like to suggest a wider range of 
possibilities: 

• Entropy history. Our simplistic schemes for 
entropy modification treat clusters as if they were 
assembled very early in time, with some standard 
initial entropy distribution. In fact, the entropy 
history of the intracluster medium is likely to be 
far more complex, influenced in differing ways at 
various times by merger shocks, radiative cooling, 
and feedback. The dispersion seen in a could stem 
from differences in the entropy history of clusters, 
with some evolutionary paths leading to abundant 
gas below the cooling threshold and other paths 
leading to very little. If this is the case, then a 
cluster's value of a would be an important window 
into its past. 

• Incomplete relaxation. Chandra observations of 
several clusters have revealed "cold fronts" thought 
to be evidence of low-entropy gas sloshing around 
in the cores of those clusters (Markevitch et al. 
2000; Vikhlinin, Markevitch, & Murray 2001). 
Our models assume that the lowest-entropy gas 
in a cluster has already settled into hydrostatic 
equilibrium at the cluster's center. Some of 

the apparent dispersion in a could be due to 
differences in the degree to which low-entropy gas 
has settled within cluster cores following the last 
major merger. 
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• Recent feedback. Cooling and condensation of 
low-entropy gas at the center of a cluster is likely 
to lead to star formation and perhaps enhanced 
activity in the nucleus of the central galaxy 
(Fabian 1994). Both star formation and radio jets 
can provide feedback, imparting entropy to the 
cluster core that inhibits further condensation 
(Bohringer et al. 2002). Perhaps these feedback 
processes control the value of a. 

• Thermal conduction. The presence of significant 
temperature gradients within the cores of some 
clusters has often been cited as evidence that 
magnetic fields strongly suppress electron thermal 
conduction in clusters (Fabian 1994). Recent 
Chandra observations of some particularly sharp 
cold fronts have reinforced that supposition (Ettori 
& Fabian 2000). However, these observations 
constrain conduction in only one dimension and do 
not necessarily rule out conduction in the other 
two dimensions. Furthermore, recent theoretical 
work suggests that magnetic fields might not 
suppress conduction as effectively as previously 
believed (Malyshkin & Kulsrud 2001; Malyshkin 
2001; Narayan & Medvedev 2001). It is possible 
that the magnetic field geometry in a cluster's core 
places a lower limit on the value of a. 

Investigating all of these possibilities is beyond the scope 
of this paper, but in the next section we explore how one 
might self-consistently calculate the intracluster entropy 
distribution within a hierarchical merging scenario. 

4. COUPLING ENTROPY, EVOLUTION, AND FEEDBACK 

We have shown that some very simple prescriptions for 
the entropy distribution in present-day clusters lead to a 
remarkably realistic set of cluster models. However, the 
astrophysics that leads to these entropy distributions is 
certainly more complicated than we have assumed. Un- 
derstanding how the the physics of structure formation, 
cooling, and feedback determine the intracluster entropy 
distribution will require a more complete theory based on 
hierarchical merging. This section outlines a few of the 
basic issues involved in constructing such a theory. 

4.1. Entropy Histories 

Considering the entropy history of a single gas parcel 
affords some insight into how radiative cooling establishes 
an entropy threshold (Voit & Bryan 2001). Suppose that 
the parcel begins with entropy K\ in a halo of tempera- 
ture T\ at some early time t\. In a hierarchical merging 
scenario, the parcel's halo will eventually merge with an- 
other halo at time ti- Merger shocks will then raise the 
parcel's entropy to some new value K 2 , and as that parcel 
settles within the new halo, it will approach a new tem- 
perature T2, similar to the characteristic temperature of 
that new halo. Each subsequent merger will also raise the 
parcel's temperature and entropy. Thus, the trajectory 
of a parcel through the entropy-temperature plane can be 
schematically represented by a sequence of points like that 
in Figure 29. 

The dashed lines in Figure 29 show how the entropy 
threshold associated with cooling rises with time until 



it achieves its present value, indicated by the solid line. 
As long as a parcel's trajectory through the entropy- 
temperature plane remains above the cooling threshold, 
then radiative cooling and subsequent feedback will not 
significantly affect that parcel's entropy. However, any 
parcel that spends a large fraction of time below that 
threshold will be subject to substantial entropy loss, con- 
densation, and whatever feedback ensues. Feedback can 
continue as long as there are gas parcels below the cooling 
threshold, but once cooling and feedback have eliminated 
all gas below the cooling threshold, feedback must cease. 
Thus, the threshold for entropy modification established 
by both cooling and feedback corresponds to the cooling 
threshold defined by the Hubble time tn ■ 

Although cooling sets the threshold for entropy modifi- 
cation, the model we are describing is not a cooling-flow 
model. Hierarchical structure formation will ensure that 
mergers disperse the products of condensation throughout 
the final cluster. Stars and galaxies form out of the lowest- 
entropy gas long before the cluster itself forms. Hence, 
much of the low-entropy gas that would eventually find 
its way into the cluster core in a simulation with no radia- 
tive cooling decouples from the intracluster medium before 
it can participate in a centrally-focused cooling flow. Fi- 
nally, for reasons we will describe next, a certain amount 
of feedback is needed to prevent the incipient intracluster 
medium from overcooling at high rcdshift. 
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Fig. 29. — Entropy-temperature diagram showing the entropy his- 
tory of a gas parcel and the evolution of the cooling threshold with 

time. Solid and dashed lines show the locus in the T-Tn e plane 
at which the cooling time t c of a gas parcel equals 15 Gyr (solid 
line), 10 Gyr (short-dashed line) and 5 Gyr (long-dashed line). This 
locus moves upward through the diagram as the universe ages. The 
dots connected by dotted lines schematically illustrate the entropy 
history of a gas parcel as described in the text. As long as the 
parcel's trajectory through the entropy-temperature plane remains 
above the cooling threshold, it will not substantially cool. However, 
both cooling and feedback triggered by that cooling will inevitably 
modify the entropy of gas parcels that find themselves below the 
cooling threshold. 



20 



MODIFIED-ENTROPY MODELS 



4.2. Overcoming Overcooling 

Overcooling is one of the classic problems that models 
of hierarchical structure formation must overcome (Cole 
1991; White & Frenk 1991; Blanchard, Valls-Gabaud, & 
Mamon 1992; Balogh et al. 2001). If radiative cooling 
were allowed to proceed unchecked as structure formed, 
then a substantial percentage (> 20%) of the universe's 
baryons would have been locked into condensed objects 
before the era of cluster formation. Because the fraction 
of uncondensed baryons in present-day clusters is <~ 90%, 
feedback processes like supernova heating are presumed to 
inhibit cooling and condensation of high-rcdshift baryons 
(White & Rees 1978; Cole 1991; White & Frenk 1991). 

In our models, overcooling manifests itself as an entropy 
threshold K c that is much higher than the characteris- 
tic entropy of an unmodified high-rcdshift halo. To illus- 
trate this effect, let us assume that the progenitor halo 
of a present-day cluster has a temperature T(z) that de- 
pends on redshift. The entropy threshold for that halo 
is then K c (z) cx {T 1 ' 2 ^) A[T(z)} t H {z)} 2 / 3 , and the char- 
acteristic entropy within r2oo in the unmodified halo is 
K 2 oo(z) °c T(z)H~ 4 / 3 . The ratio of these two quantities 
in the high-redshift limit is thus 

P~ « l^TT 1 2 ' 3 mz)t H (z)f/ 3 H^(z) . (22) 
if 200 I T(z) J 

Because H(z) rises and T(z) declines with increasing red- 
shift, this ratio must exceed unity at early times. 

We can use extended Press-Schechter theory (Bond et 
al. 1991; Bower 1991) to approximate T(z) and esti- 
mate the temperature, entropy, and redshift scales at 
which overcooling is problematic. First, to demonstrate 
the principle with simple scaling relations, we consider a 
SCDM (fijkf = 1) universe with a power spectrum given 
by an n = —2 power law. In that case, the characteris- 
tic mass scale evolves like M cx (1 + z)~ 6 , and assuming 
that M cx T 3 / 2 (l + z)- 3 ' 2 yields T(z) cx (1 + z)~ 3 . At 
low redshifts, where T is high, we have A cx T 1 / 2 and 
K c /K 2 oo cx (1 + z) 2 . However, at high redshifts, when 
the progenitor temperature is smaller, we have A cx T~ x l 2 
and K c /K 200 cx (1 + z) 4 . Thus, the progenitor of a halo 
with K c /K 2no ~ 0.1 at z m would have K c /K 200 ~ 1 at 
z ~ 1 - 2. 

Figure 30 illustrates this effect more precisely for a 
ACDM universe seeded by a power spectrum with shape 
parameter T = 0.17. The dashed line shows K 2 oq{z) for 
a halo that reaches T — 5 keV at z = 0. The solid line 
shows the entropy threshold K c in that evolving halo for 
a metallicity of one-third solar. Feedback must impart an 
entropy ~ 30keVcm 2 to prevent overcooling at z > 3. 

4.3. Interplay between Feedback and Merging 

Feedback can solve the overcooling problem, but it in- 
troduces another, more subtle problem into our frame- 
work. Our modified-entropy models assume that gas above 
the cooling threshold follows the unmodified distribution 
K (M g ). However, the need for feedback at early times 
implies that a large proportion of intracluster gas has ex- 
perienced entropy modification at least once in its history. 
Somehow, the process of cluster formation needs to restore 
the unmodified distribution to gas above the threshold at 
z = 0. 




(1 + z) 

Fig. 30. — Entropy evolution with redshift. The dashed line 
shows how the characteristic entropy K200 of the progenitor halo 
of a 5keV cluster at 2 = evolves with redshift in a ACDM uni- 
verse with shape parameter T = 0.17. The solid line shows how the 
cooling threshold K c associated with that halo depends on redshift. 
Because K200 falls below the cooling threshold at z > 3, feedback is 
needed to prevent the majority of the baryons associated with the 
progenitor from condensing and forming stars. 



Exactly how that happens is unclear. Consider what 
happens in a merger between two halos containing isother- 
mal gas at temperatures T\ and T 2 that collide at velocity 

— 1/2 

v. Merger shocks with Mach number Mi cx vT i will 
propagate through each halo's gas, raising the entropy of 
that gas by some more-or-less uniform factor depending on 
Mi. This mechanism could potentially amplify the effects 
of early feedback, which must raise the intrahalo entropy 
to at least K c to avoid overcooling. Subsequent shocks 
would then boost the entropy even further, along a path 
parallel to the dashed line in Figure 30. If this mechanism 
dominates, then very little gas should remain near K c at 
z~ 0. 

However, early feedback is likely to diminish the effects 
of mergers, particularly in gas with low entropy. If feed- 
back raises all the gas in the progenitor halo of a clus- 
ter to ~ K c , then the intrahalo medium will be nearly 
isentropic, with high-temperature gas at the center of the 
halo and lower-temperature gas at the outskirts (see Fig- 
ure 8). When mergers occur, Mach numbers associated 
with lower-entropy gas in the halo's core will therefore be 
lower than those associated with higher-entropy gas near 
the virial radius. Hence, merger shocks will act to steepen 
any entropy gradient that exists. Dynamical friction may 
also add to this effect, because the cores of merging halos 
tend to lose orbital energy before they collide and merge. 
Another effect to consider is the quasi-continuous accre- 
tion of gas associated with small halos, for which the Mach 
number is likely to be quite high. This type of accretion 
will enhance the amount of high-entropy gas in the final 
cluster. Thus, it seems possible that the modified-entropy 
models we have developed can successfully be linked with 
hierarchical structure formation, although many details re- 
main to be worked out. 
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4.4. Excess Energy 

Several recent analyses of similarity breaking in groups 
and clusters have cast the problem in terms of energy 
rather than entropy (e.g., Wu et al. 1998, 2000; Lloyd- 
Davies, Ponman, & Cannon 2000; Loewenstein 2000; 
Bower et al. 2001; Lloyd-Davies et al. 2002). For ex- 
ample, if one knows what the configuration of the intra- 
cluster gas would be in the absence of non-gravitational 
heating and cooling processes, then one can opt to define 
the "excess energy" of a cluster's gas to be the difference 
in total energy between that unmodified configuration and 
the actual configuration. Depending on what is assumed 
about the unmodified configuration, these excess energies 
can range from ~ 0.3 — 1.0 keV for groups (e.g., Lloyd- 
Davies et al. 2000, 2002) to over 2 keV for clusters (e.g., 
Wu et al. 2000). 

The results of this paper indicate that some combination 
of heating and cooling is responsible for similarity break- 
ing because cooling is needed to explain the core entropies 
of clusters and groups and heating is needed to prevent 
overcooling at early times. Therefore, we prefer to com- 
pute excess energy in a way that explicitly accounts for the 
separate contributions of cooling and heating. First, we as- 
sume that the baseline state of a cluster of mass M200 is a 
truncated model from which cooling and condensation has 
removed a fraction /* of the lowest-entropy gas. In other 
words, the initial configuration of the intracluster medium 
has an entropy floor equal to K (f*) and contains a frac- 
tion 1 — /* of the cluster's baryons. Then, we assume that 
the final state of the cluster is a shifted model with the 
same proportion of baryons but an entropy floor equal to 
the cooling threshold, so that it obeys the observed M-T 
and L-T relations. Formally, the entropy distribution of 
the final state is K(f g ) = Ko(f g ) + K c (T\ um ) with an outer 
boundary where f g = 1 — /* . Because the amount of hot 
baryonic gas in both models is the same, the amount of 
heat input required to explain the cluster scaling relations 
after a fraction /* of the intracluster gas has condensed 
equals the energy input needed to convert the initial con- 
figuration to the final configuration, including work done 
at the outer boundary of the intracluster medium. 

Figure 31 shows the energy difference between these con- 
figurations as a function of halo mass for different values 
of the condensed baryon fraction /*. To facilitate com- 
parisons with other treatments of excess energy, we di- 
vide this energy difference by the total number of parti- 
cles in an unmodified cluster (/hM2oo/M m p)j gi ym g the 
quantity Ae in units of keV per particle. When the con- 
densed baryon fraction is small (/» = 0.02), the excess en- 
ergy needed to account for the scaling relations of massive 
clusters (~ 10 15 M Q ) approaches the ~ 3keV level. How- 
ever, this large amount of excess energy depends heavily 
on the assumed form of the unmodified entropy distribu- 
tion, which is likely to be a poor approximation for these 
purposes. The large difference in the excess energy of mas- 
sive clusters between models with /* = 0.02 and /* > 0.05 
indicates that a large amount of energy is needed to lift the 
lowest-entropy gas out of the very center of the gravita- 
tional potential. Yet, Figure 1 shows that the unmodified 
intracluster density distribution is unlikely to be as cen- 
trally condensed as we have assumed. Even in the absence 
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Fig. 31. — Excess energy required to overcome the cooling thresh- 
old for different values of the condensed baryon fraction (/*). At 
each halo mass M20O! the quantity Ae[-fsT c (Ti um ), /»] is the energy in- 
put per particle needed to transform the truncated entropy distribu- 
tion Kt (/* ) , from which a fraction /* of the lowest-entropy baryons 
have been removed, into a shifted distribution Kg = Kq + K c 
containing the remaining fraction 1 — /* of the baryons. This 
transformation corresponds to first removing gas without adding 
heat energy, then adding enough heat to reproduce the L-T and 
M-T relations. The energy input required to produce these rela- 
tions for plausible condensed-baryon fractions of /„ fs 0.1 — 0.2 is 
~ 0.4 — 0.7 keV particle -1 at group scales and is negligible for hot 
clusters. 

of cooling, the true unmodified profile probably corre- 
sponds more closely to a truncated model with K c ~ 
0.05 — 0.1, which substantially reduces the excess energy 
burden. Thus, when calculating excess energies for clus- 
ters, one must make sure that the baseline model accu- 
rately represents the cluster configuration in the absence 
of heating. 

The curve with /* = 0.1 is likely to be a more realis- 
tic representation of the heat input needed to explain the 
scaling properties of massive clusters for two reasons: (1) 
this value of /* roughly corresponds to the fraction of clus- 
ter baryons that have condensed into stars, and (2) when 
/* ^ 0-1) truncation eliminates the innermost part of the 
entropy distribution that is discrepant with simulations. 
At 10 15 Mq, this excess energy curve suggests that only a 
fraction of a keV per particle is needed to explain cluster 
scaling relations. That is because a shifted model with an 
entropy shift determined by the cooling threshold is virtu- 
ally identical to a truncated model with the innermost /* 
of gas removed. Likewise, the mass scales at which curves 
with larger values of /* intercept the zero-energy level also 
indicate cases in which the truncated and shifted models 
are virtually identical, requiring no excess heating to pro- 
duce the observed scaling relations. In the limiting case 
of zero supernova heating (e.g., Bryan 2000, Muanwong et 
al. 2001), these intercept points show the fraction of con- 
densed baryons implied by the L-T relation as a function 
of halo mass. 

If the condensed baryon fractions of groups are similar 
to those of clusters, then <~ 0.3 — 0.7 keV of excess energy is 
needed to explain the M-T and L-T relations at ~ 1 keV. 
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The main difference between the initial configuration and 
the final configuration of groups with /* ~ 0.1 — 0.25 is 
that the gas in the final configuration is considerably more 
extended, with substantially more gravitational potential 
energy. Because models with higher /» contain less in- 
tragroup gas, they are less extended, and less energy is 
needed to change their configuration. We note that the 
excess energy needs of groups are comparable to estimates 
of the supernova energy released over the course of the 
group's history (e.g., Ponman et al. 1999, Loewcnstcin 
2000; Pipino et al. 2002), which may mean that no other 
source of energy is needed to explain their scaling proper- 
ties, if that energy efficiently heats the intragroup medium. 

The bottom line here is that calculations of the heat 
input needed to account for the scaling relations of clus- 
ters and groups depend critically on what is assumed to 
happen in the absence of heat input. Most calculations 
of excess energy have assumed that the intracluster gas 
density obeys either a simple polytropic equation of state 
(e.g., Wu et al. 2000; Loewenstein 2000) or a /3-model 
density profile (e.g., Bower et al. 2001; Lloyd-Davies et 
al. 2002) whose parameters are adjusted according to the 
amount of heat input. The modificd-entropy models dis- 
cussed in this paper allow for cooling and condensation of 
the lowest-entropy gas to occur before heat is added. Fig- 
ure 31 demonstrates that the amount of heating needed 
to explain the scaling relations depends strongly on the 
amount of cooling that occurs prior to heat input and that 
realistic amounts of "precooling" substantially reduce the 
required amount of non-gravitational heating. 

5. CONCLUSIONS 

A realistic family of models for clusters in hydrostatic 
and convective equilibrium can be constructed using the 
Navarro, Frenk, & White (1997) density profile and some 
rudimentary prescriptions for how radiative cooling and 
feedback induced by that cooling modify the entropy dis- 
tribution of intracluster gas. The models presume that the 
intracluster entropy distribution in the absence of cooling 
and feedback would be identical to that of gas that takes 
on the same density distribution as the dark matter as 
it settles into the cluster's potential well — an assumption 
supported by numerical simulations. Because the lowest- 
entropy gas within that distribution can cool within a Hub- 
ble time, its entropy must somehow be modified. The 
threshold entropy in our prescriptions for entropy modi- 
fication corresponds to the entropy at which the cooling 
time of intracluster gas equals the age of the universe. The 
prescriptions themselves include truncation of the entropy 
distribution at the threshold entropy, shifting of the en- 
tropy distribution by adding the threshold entropy to the 
entire distribution, and a qualitative implementation of 
radiative cooling. 

Exploring the properties of dimensionlcss models de- 
pending only on the halo concentration c, the dimension- 
less entropy threshold K c , and the prescription for entropy 
modification reveals that: 

• Removal of low-entropy gas acts to flatten the 
core density profile, regardless of whether cooling 
or heating eliminates that low-entropy gas (see 
Figures 1 and 4). 



• The observable properties of clusters depend 
more critically on the threshold entropy than on 
the mode of entropy modification. For example, 
truncation of the intracluster entropy distribution 
at a given value of K c and shifting of the entropy 
distribution by that same amount lead to very 
similar pressure, density, and temperature profiles 
within the cluster's virial radius (see Figures 4 and 
8). 

• Within about 30% of the virial radius, the 
surface-brightness profiles of truncated and shifted 
models are very similar to /3-models (see Figure 5). 
Raising the entropy threshold tends to lower the 
best-fitting value of 0, while raising the halo 
concentration tends to decrease the core radius 
(see Figures 6 and 7). 

• All of our prescriptions for entropy modification 
tend to augment the luminosity-weighted temper- 
ature of a cluster as the entropy threshold rises. 
However, that temperature increase is relatively 
small because the increased central entropy 
decreases the luminosity of high-temperature gas 
in the core. Thus, modifiying the entropy of a 
cluster shifts the bulk of a cluster's luminosity 
to larger radii, where gas temperatures tend to 
be smaller, mitigating the effects of the entropy 
increase (see § 2.5). 

• All of our prescriptions for entropy modification 
reduce the X-ray luminosity of a cluster as the 
entropy threshold rises. In the high-threshold limit, 

dimensionless luminosity scales as L cx i\ f T 3/ ' 2 , 
owing to the asymptotic p cx r~ 3 scaling of the 
unmodified NFW density distribution at large 
radius (see § 2.6). When luminosity scales with 
the entropy threshold in this way, the L-T relation 
should not flatten below ~ 2keV, where line 
cooling dominates free-free emission, and the 
relation should also evolve little with time. The 
observed L-T relation shares both of these features. 
Furthermore, this scaling of luminosity with the 
entropy threshold leads to an L-T relation slightly 
steeper than L cx T lu ' m , with the extra steepening 
coming from the tendency for low-temperature 
halos to be more concentrated, also in agreement 
with observations. 

Adopting relations between halo concentration c and 
halo mass M200 that are drawn from simulations and sup- 
ported by observations enables us to generate models for 
real clusters that reproduce many of their observable prop- 
erties: 

• Fitting /3-models to our modified-entropy clusters 
yields (3 and core-radius values that are similar to 
those of observed clusters. Our models reproduce 
the observed tendency for (3 to be lower in low- 
temperature halos, a behavior that arises because 
the entropy threshold determined by cooling has a 
much larger impact on low-temperature halos. Our 
models also reproduce the observed relationship 
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between core radius and temperature, implying 
that r c w 0.1r2oo, regardless of temperature. 
Despite the flattening of the density profile in 
low-temperature halos, the core radius remains 
near 0.1r2oo because of the increased halo 
concentration (see Figures 6 and 7). The large 
dispersion observed in (3 and r c could arise from 
a number of effects, some involving observational 
systematics such as surface-brightness bias and the 
range of radii in the fit, others involving physical 
differences such as cluster-to-cluster variations in 
the amount of gas below the cooling threshold (see 
Figures 17 and 18). 

• The temperature gradient at small radii in our 
model clusters depends sensitively on the amount 
of gas below the cooling threshold, closely related 
to the parameter a in our radiative-loss models. 
In our most extreme radiative-loss model, with 

a = 3/2, the temperature gradient at < 0.3 r 200 
is very similar to the "universal" temperature 
gradient observed by Allen et al. (2001). As the 
parameter a rises to infinity, the limit in which 
no gas lies below the cooling threshold, the inner 
temperature gradient flattens and then reverses, 
monotonically decreasing with radius in the large 
a limit (see Figure 19). 

• All of our models have a negative temperature 
gradient at large radii because the underlying 
potential is steeper than isothermal at those 
radii. In order to compare our models with 
observations, we compute the effective polytropic 
index 7 e ff = d In P/d In p as a function of radius. At 
the virial radius, we find 7 e g w 1.1 — 1.2, but the 
behavior of 7 e g at smaller radii depends strongly 
on the nature of entropy modification. In models 
with no gas below the cooling threshold, the 
cluster core is nearly isentropic, with 7 e g- w 5/3. 
In models with significant amounts of gas below 
that threshold, we find 7 c g- < 1 within the core 
radius, a consequence of the positive temperature 
gradient (see Figure 20). Values of 7 e ff evaluated 
at 0.2 r2oo agree with observations, within the large 
observational uncertainties (see Figure 21). 

• The mass-temperature relations derived from 
our truncated and shifted modificd-entropy 
models, which have no gas below the cooling 
threshold, agree well with those derived from 
cluster observations that have been corrected for 
the presence of cooling flows (see Figure 22). The 
slope of this relation is steeper than the M oc T 3 / 2 
expectation from self-similar scaling because 
low-temperature halos are more concentrated, 
leading to a slightly higher temperature for a given 
mass (see Figure 9), a feature these models share 
with those of Komatsu & Seljak (2001). Clusters 
of a given temperature are less massive than 
those in simulations without cooling and feedback 
because entropy modification owing to these 
processes shifts all of our model clusters to higher 
temperature. However, the agreement between 



observations and our models does not necessarily 
imply that this M-T relation is without problems, 
as both the observational determinations and the 
models assume that the intracluster medium is an 
ideal gas in hydrostatic equilibrium, which may 
not be the case. Furthermore, this relation does 
not apply to clusters that have not been corrected 
for cooling flows — the temperature normalization 
at a given mass in our radiative-loss models can be 
up to 30% lower, depending on the value of a (see 
Figures 11 and 12). 

• The luminosity-temperature relations derived 
from truncated and shifted models also agree well 
with those derived from cooling-flow corrected 
observations of clusters (see Figure 25). While 
there are no adjustable parameters in these models, 
the baryon-to-dark-matter ratio does depend on 
our assumed value of Q,m = 0.33. Thus, this 
agreement can be taken as additional evidence for 
a low value of Om- 

• Gas below the cooling threshold both raises the X- 
ray luminosity and lowers the luminosity-weighted 
temperature of our model clusters. Hence, the L-T 
relations derived from our radiative-loss models lie 
above and to the left of the truncated and shifted 
models in the L-T plane. Changing a leads to 

a displacement almost perpendicular to the L-T 
relations themselves (see Figure 28). Correcting 
for cooling-flow emission leads to a very similar 
displacement in L-T space, suggesting that the 
L-T relation for uncorrected clusters may depend 
almost entirely on the two parameters M 2 oo and 
a. It will be interesting to see whether values of a 
derived from displacement in L-T space agree with 
those implied by the inner temperature gradient. 

The modified-entropy models we have derived can ac- 
count for many of the observed properties of present-day 
clusters, but they are obviously too simplistic to tell us 
how those clusters got to be the way they are. Somehow, 
hierarchical structure formation generated intracluster en- 
tropy distributions similar to our simple prescriptions for 
entropy modification. Cooling, feedback, and merging all 
have important roles to play. The consistency of our mod- 
els with the data strongly suggests that radiative cooling 
sets the level of the intracluster entropy threshold and sub- 
stantially reduces the amount of heat input needed to ex- 
plain the observed scaling relations. However, radiative 
cooling unchecked by feedback would lead to overcooling 
of intrahalo gas at z <~ 2 — 3 (see § 4.2). At least some feed- 
back is needed to keep ~ 90% of intergalactic baryons in 
gaseous form, as observed. Mergers following an episode of 
feedback can potentially amplify the entropy generated by 
that feedback, when low-density gas associated with low- 
temperature halos is shocked. However, the overall impact 
of mergers on intracluster entropy remains unclear. Deter- 
mining just how cooling, mergers, and feedback conspire 
to produce the present-day entropy distribution of clus- 
ters will require high-resolution numerical simulations and 
semi-analytical modeling focusing on the intergalactic en- 
tropy distribution and how it evolves. 
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